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ABSTRACT 

This  work  addresses  the  design  and  analysis  of  the 
Pitch,  Yaw  and  Roll  autopilot  for  application  to  the 
Bank-To-Turn  (BTT)  missiles.   At  first,  the  linear 
uncoupled  channels  were  designed  and  analyzed  according  to 
the  desired  requirements.   Utilizing  the  uncoupled 
channels,  the  linear  coupled  autopilots  were  designed,  not 
including  the  inertial,  kinematic  and  aerodynamic 
cross-coupling.   Then,  the  nonlinear  CBTT  autopilots  were 
designed  and  analyzed,  using  the  linear  CBTT  (coordinated 
Bank-To-Turn)  models,  which  now  have  coupled  with 
kinematic,  inertial  and  aerodynamic  cross-coupling.   The 
minimization  of  the  above  kinematic  and  inertial  coupling 
and  their  effects  were  completed  using  feedbacks  of 
angle-of-attack  and  rate  of  angle-of-attack  in  the  Pitch 
autopilot. 
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TABLES  OF  SYMBOLS  AND  ABBREVIATION 


BTT 
CBTT 

CRDT 
C. 


Bank-to-Turn 

coordinated  Bank-to-Turn,  minimum  sideslip, 
positive  a,  4     180  deg 

coordinated  Roll-During-Turn 
rolling  moment  coefficient 

slope  of  curve  of  rolling  moment  coefficient, 

C   vs  B 

I 


change  in  C   per  degree  roll  control  incidence, 


change  in  C   per  degree  yaw  control  incidence, 
5Y 


m 


m. 


m 


'N 


'N 


pitching  moment  coefficient 

slope  of  curve  of  pitching  moment  coefficient 

C   vs  a 
m 

change  in  C   per  degree  pitch  control 

incidence,  6_ 

normal  force  coefficient 

slope  of  curve  of  normal  force  coefficient  C 


vs  a 


N 


'N 


change  in  CN  per  degree  pitch  control 
incidence,  6p 


n 


'n 


yawing  moment  coefficient 

slope  of  curve  of  yawing  moment  coefficient,  Cn 
vs  b 
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q  pitch  angular  acceleration  about  YB 

Q  constant  or  equilibrium  pitch  angular  rate 


r  yaw  angular  rate  about  Z 


B 


rc  yaw  angular  rate  command  (coordination  command) 


B 

2 


r  yaw  angular  acceleration  about  Z 

S  reference  area  for  coefficients  =  ir  ft 

STT         Skid-to-Turn,  roll  attitude  stabilized 

u  velocity  component  in  XB  direction 

v  velocity  component  in  Y_  direction,  assumed  to 

be  constant 

V  constant  missile  flight  path  velocity 

V  missile  velocity  vector 

w  velocity  component  in  ZR  direction 

X„  body-fixed  roll  axis,  along  axis  of  symmetry, 

positive  forward 

V  body-fixed  pitch  axis,  positive  starboard 

yv  vehicle  axis  in  local  horizontal  direction, 

approximated  as  inertial  axis 

zB  body-fixed  yaw  axis,  forms  right  handed 

orthogonal  system  with  XB  and  Yg 

z  vehicle  axis  in  downward  direction  along  local 

gravity  vector,  approximated  as  inertial  axis 

n  achieved  normal  acceleration  in  Z„  direction 

z  V 

nz  commanded  normal  acceleration  in  ZB  direction 

c 

n  achieved  normal  acceleration  in  YR  direction 

yB 

n„  achieved  normal  acceleration  in  ZD  direction 

ZB  B 
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C  change  in  C   per  deqree  yaw  control  incidence. 

6Y  6Y 

C  change  in  C   per  degree  yaw  control  incidence, 

^R  SR 

C  side  force  coefficient 

C  slope  of  curve  of  side  force  coefficient  C   vs 

6  8 

CY  change  in  Cy  Per  degree  yaw  control  incidence, 

CY  change  in  CY  Pe*"  degree  roll  control  incidence, 

6R  6R 

d  reference  length  for  cefficients  =  2  ft. 

I  moment  of  inertia  about  Y_  axis 

yy  B 

Izz  moment  of  inertia  about  ZB  axis 

I  moment  of  inertia  about  X_  axis 

XX  b 

K^  autopilot  pitch  acceleration  error  gain 

Kyp  CBTT  autopilot  coordination  branch  gain 

LBTT  limited  Bank-to-Turn  may  or  may  not  be 

coordinated,  positive  and  negative  a,  <t>  ^  90 
or  45  degrees 

P  roll  rate  about  X 

P  roll  acceleration  about  XB 

P  constant  or  equilbrium  roll  angular  rate 

POC  preferred  orientation  control 

q  dynamic  pressure 

q  pitch  rate  about  YB 
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r\y  achieved  normal  acceleration  in  Yv  direction 

n  normal  acceleration  command  from  guidance 

computer  in  Zv  direction  plus  anti-gravity  bias 
command 

<fr  roll  attitude  command  from  guidance  computer, 

zero  degrees  in  Z..  direction  and  90  degrees  in 
Y   direction 

<j>  roll  attitude,  zero  degrees  in  -Zv  direction 

and  90  degrees  in  Yv  direction 

e  Elevation  Euler  Angle,  second  rotation, 

(q  cos<|>  -  r  sin4»)  dt 

5_  pitch  control  incidence  (positive  tail 

incidence  produces  negative  pitching  moment)* 

5  commanded  pitch  control  incidence,  5p 

c 

«Y  yaw  control  incidence  (positive  tail  incidence 

produces  negative  yawing  moments)* 

6y  commanded  yaw  control  incidence,  5„ 

c 

5  roll  conntrol  incidence  (positive  tail 

incidence  produces  positive  rolling  moment) 

6„  commanded  roll  control  incidence,  6R 

c 

a  constant  or  equilibrium  angle-of-attack 

a  angle-of-attack 

a  ange-of-attack  rate 

a  modified  form  of  estimated  angle-of-attack  for 

autopilot  coordination  command 

8  angle  of  sideslip 

< 

B  sideslip  angular  rate 
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I.   INTRODUCTION 

In  recent  years,  the  application  of  Bank-to-turn 
guidance  to  missiles  has  generated  considerable  interest. 
This  interest  has  been  motivated  by  the  drag  reduction, 
which  the  Bank-to-turn  (BTT)  controlled  missiles  offer, 
over  the  conventional  cruciform  roll  stabilized, 
skid-to-turn  controlled  missiles. 

Bank-to-turn  steering  may  provide  improved  performance 
for  missile  systems.   In  addition,  the  ramjet  engine  inlet 
configuration  is  believed  to  provide  a  greater  range 
capability  than  other  inlet  configurations.   BTT  control  is 
required  to  satisfy  the  sideslip  constraints  imposed  by 
chin  inlets.   Despite  these  facts,  there  are  unanswered 
questions  concerning  BTT  systems  stability  during  homing, 
guidance  performance,  autopilot  and  guidance  logic  design 
and  subsystem  requirements  which  must  be  investigated 
before  Bank-to-turn  steering  can  be  considered  a  viable 
method  of  control  for  high  performance  missiles. 

Many  missile  programs  were  initiated  during  the  past 
decade  to  improve  their  capability  via  Bank-to-turn 
control.   The  results  have  advanced  the  understanding  of 
the  different  missile  subsystems.   In  the  autopilot  area, 
many  types  of  autopilots  have  been  found  which  force  the 
missile  to  roll  or  bank  so  that  the  steering  maneuver 
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occurs  with  the  airframe  oriented  in  a  specified  or 
preferred  direction  with  respect  to  the  incoming  airstream. 
This  entire  class  of  autopilots  may  be  referred  to  as 
preferred  orientation  control  autopilots. 

Also,  simplified  guidance  studies  which  neglect  radome 
effects  and  assume  coordinated  missile  motion  have  shown 
that  CBTT  can  provide  acceptable  performance  with  roll 
rates  that  are  not  excessive  for  autopilot  control.   These 
studies  were  made  for  a  medium  range  area  defense  mission 
and  a  long  range  suppression  mission  and  considered  both 
high  lift  and  moderate  lift  configurations  [Ref.  1] . 

The  present  work  reviews  the  design  and  analysis  of 
the  pitch,  Yaw  and  Roll  autopilot  for  application  to  the 
Bank-to-turn  (BTT)  missiles  developed  in  [Ref,  2].   At 
first,  the  linear  model  of  the  uncoupled  autopilot  channels 
will  be  designed  and  analysed  in  the  critical  areas  of 
concern  regarding  the  response  of  CBTT  control  for  the 
circular  and  elliptical  airframes.   Then,  the  designed 
autopilot  channels  will  be  coupled  using  suitable 
aerodynamic  models  and  a  measure  of  sideslip  control  is 
obtained  by  roll  command  to  the  linearized  CBTT  autopilot. 
The  relative  stability  of  the  autopilot  branches  and  means 
for  improving  stability  will  be  discussed.   In  a  later 
chapter  the  non-linear  three  dimensional  model  will  be 
designed,  using  a  three  dimensional  aerodynamic  model  and 
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the  same  control  laws  described  in   the  linear  studies, 
except  for  a  minor  modification  to  the  coordinating  branch 
dependence  on  Angle-of-Attack  and  also  the  inclusion  of 
Anti-Gravity  Bias.   Next,  the  kinematic  and  inertial 
coupling  effects  between  pitch  and  yaw  dynamics  will  be 
studied.   Transients  in  maneuver  plane  acceleration  are 
caused  in  the  above  cross-coupling  effects  and  are  in  the 
form  of  overshoots  and  undershoots.   Coupling  transients 
will  be  reduced  by  increasing  pitch  stability.   The 
technique  of  feedbacks  of  Angle-of-Attack  and  rate  of 
Angle-of-Attack  for  the  pitch  autopilot,  will  be 
considered,  in  order  to  decrease  kinematic  and  inertial 
cross-coupling  effects.   Finally  conclusions  and 
recommendations  for  future  study  will  be  stated. 
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II.   LINEAR  DESIGN  AND  ANALYSIS  OF  UNCOUPLED  AUTOPILOTS 
FOR  CIRCULAR  AND  ELLIPTICAL  AIRFRAME 

A.    GENERAL 

The  linear  design  and  analysis  started  with  uncoupled 
autopilot  channels,  i.e.  aerodynamic,  kinematic  and  control 
law  cross-coupling  between  Pitch,  Yaw  and  Roll  channels 
were  removed.   The  autopilot  design  technique  was 
classical,  using  a  combination  of  frequency  response  and 
root  locus  techniques  [Ref .  2] ,  providing  the  range  of 
required  missile  body  angular  rates,  control  motions  and 
minimizing  the  influence  of  aerodynamic  variations  on 
desired  responses.   To  achieve  the  desired  maneuver  plane 
acceleration  response  for  the  CBTT  autopilot,  both  Pitch 
and  Roll  uncoupled  autopilots  were  designed  to  have  a  time 
constant  of  0.5  sec  in  order  to  achieve  the  desired 
maneuver  plane  speed  of  response.   The  Yaw  channel  which 
follows  the  Roll  channel  motion  to  produce  desired 
coordination  (i.e.  minimization  of  sideslip  angles),  was 
designed  to  have  a  more  rapid  response  than  the  Roll 
channel,  having  0.39  sec  for  circular  and  0.36  sec  for 
elliptical  time  constant  in  the  maneuver  plane  acceleration 
response.   The  whole  design  was  done  for  the  constant 
flight  condition  of  60  KFT  altitude  and  mach  number  3.95. 
The  design  discussed  here  was  developed  by  Arrow  [Ref.  21. 
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In  this  chapter  the  airframe  configurations,  the 
aerodynamic  models  and  control  laws  for  the  Pitch,  Yaw  and 
Roll  channel  are  reviewed.   The  transfer  functions  for  both 
circular  and  elliptical  configuration  are  derived.   An 
analysis  for  every  channel  in  regard  to  the  acceleration 
responses,  the  body  angular  rates  and  control  surface 
deflection  is  made.   Also,  a  comparison  of  airframe 
responses  for  every  channel  is  made,  considering  the 
relative  stability  and  required  control  surfaces. 

A  general  block  diagram,  of  the  channels  and  the 
command,  which  are  applied  to  the  control  laws  is  shown  in 
Figure  2.1.   A  number  of  the  Figures  shown  in  this  report 
have  been  taken  from  [Ref.  2]. 

B.    AIRFRAME  CONFIGURATIONS 

The  U.S.  National  Aeronautics  and  Space  Administration 
maintains  a  continuing  research  effort  in  missile  related 
aerodynamics.   In  this  effort  monoplanar  missile  concepts 
have  been  considered,  in  part,  because  of  the  relatively 
low  geometric  profile  that  can  be  achieved  for  convenient 
carriage  [Ref.  1].   The  two  Airframe  Configurations  studied 
in  this  research  were  taken  from  [Ref.  1]  and  are  shown  in 
Figure  2.2  and  Figure  2.3.   The  circular  cross  sectional 

body  has  a  closure  ratio  A.  ,/A,,,^  of  0.69  and  the  A  av 

DaSc   max  max 

occurs  at  68.0%  body.   The  Elliptical  has  an  exact  3:1 
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elliptical  cross  section,  and  its  cross  section  area  is  the 
same  as  the  area  of  circular  body.   The  span  of  both 
airframes  is  the  same.   The  wing  area  and  span  for  the 
circular  concept  were  chosen  as  typical  of  current 
maneuvering  missiles.   The  wing  for  the  elliptical 
configuration  was  determined  by  projecting  the  elliptical 
body  on  the  circular  bodywing  planform.   The  resultant 
exposed  wing  planform  then  became  the  wing  for  the 
elliptical  body.   The  tails  for  each  design  had  a  dihedral 
of   30   degrees  as  shown  in  Figure  2.2.   In  order  for  tail 
deflection  to  be  compatible  with  complex  surfaces  of  the 
afterbody  of  the  elliptical  configuration,  the  tail  hinge 
line  was  skewed  such  that  a  10  degrees  deflection  measured 
at  the  body-tail  juncture  had  a  resultant  7.04  degrees 
surface  deflection,  as  shown  in  Figure  2.3  and  [Ref.  1]. 
Comparison  of  the  elliptical  cross  section  concept  with  the 
circular  cross  section  concept  indicates  the  following: 

1.  About  25  percent  more  normal  force,  that  is  nearly 
independent  of  angle  of  attack  can  be  achieved  at 
supersonic  speeds. 

2.  Values  of  the  longitudinal  stability  parameter  C, 
are  more  positive,  with  more  pronounced  nonlinearities  in 
pitching  moment  occur ing  at  subsonic  speeds. 

3.  Levels  of  directional  stability  are  increased  and 
are  more  compatible  with  levels  of  longitudinal  stability. 
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4.   Noticeably  more  yaw  control  is  available,  although 
suitable  locations  for  tails  on  the  body  are  more  limited 
because  of  the  geometry  of  the  elliptical  body. 

C.    UNCOUPLED  LINEAR  PITCH  CHANNEL 
1.   Aerodynamic  Model 

The  initial  phase  in  the  design  of  the  CBTT 
autopilots  involves  the  design  of  the  uncoupled  pitch 
channel.   The  linear  pitch  channel  dynamic  model  which  is 
used,  is  shown  in  Figure  2.5. 

This  figure  is  taken  from  [Ref.  2].   The  above 
model  is  linearized  for  stability  studies.   The  following 
three  assumptions  were  made: 

a.  Plane  Xg-Zg  is  the  maneuver  plane,  shown  in 
figure  2.4 

b.  Missile  is  trimmed  in  pitch,  that  is  My=0,  at 
fixed  values  of  a,  q  and  5  . 

c.  Missile  roll  rate  is  constant.   Using  the 
foregoing  assumptions,  and  the  approximation  that  the 
variation  in  forward  velocity  (u),  is  equal  to  zero,  a 
simplified  set  of  longitudinal  equations  of  motion  are 
constructed  applicable  to  short  period  motion.   These 
equations  describe  the  two  degree  of  freedom,  short  period 
mode  longitudinal  motion. 

Linearized  aerodynamic  derivatives  for  M=3.95 
are  provided  in  Table  I. 


n 

M— 

3 

_a> 

0) 

n> 

to 

to 

O 

O 

c 

c 

CD 

+ 

+ 

+ 

o 

c 

03 

n 

1 

+ 

1 

■o 

«3 

1 

_ „ 

CN 

+ 

+ 

1 

ca 

- 

1 

+ 

+ 

"O 

"O 

c 

.    ro 

■<5  E 

u 

(J 

H  E 

a. 

> 

a: 

o 

>-o 

<o 

<<3 

o 

-r 

+ 

+ 

c 
o 

•cH 

4J 

o 
a> 
> 

c 
o 
u 

c 

•tH 

U 


s    • 

<u  t/i 

C   --H 

>  X 

T3    <0 

o 

U  T3 

d)  c 

<  < 

c 

•^ 

1-^- 

o 

• 

o 

*-> 

CN 

c 
o 

o 

♦^ 

CD 

CT> 

u 

a> 

TJ 

•■-t 

k_ 

CD 

Du 

13 

TJ 

«/> 

CD 

CD 

O) 

ro 

c 

C> 

T3 

U 

ca 

c 

CD 

5 

CD 

> 

o 

t_ 

*-< 

< 

en 

o 

32 


U 
O 


a. 
o 

3 
< 


(!) 

C 
C 
<0 

-C 

o 

u 

04 

<U 

rH 

a 

3 
O 

o 

c 

=3 


CM 

« 

&4 


33 


TABLE  I 


Linearized  Aerodynamic  Derivatives  (M  =  3.95) 


Circular 


Elliptical 


a  =  0°    a=  10°   a  =  20°   a  =  0°    a    =  10°       a  =  20° 


-  .065    -  .082    -  .111    -  .043    -  .054 


-  .064 


n 


-  .025    -  .019    -  .003    +  .024    +  .024 


+  .032 


-  .009    -  .020 


-  .027 


-  .040 


'Y 


n 


+  .021    +  .022    +  .028    +  .016    +  .015 


-  .050    -  .053    -  .062    -  .042    -  .039 


+  .019 


-  .045 


-  .016    -  .038 


-  .010 


-  .023 


R 


n 


-  .009    -  .022 


+  .018    +  .044 


-  .006 


+  .014 


-  .014 
+  .032 


'R 


'N 


+  .031    +  .035    +  .044    +  .023    +  .023 


+  .15     +  .17     +  .22     +  .18     +  .22 


+  .029 


+  .30 


'N 


+  .04     +  .04     +  .05     +  .02     +  .02 


+  .025 


m 


-  .060    -  .065    -  .118    +  .015    +  .0137      -  .125 


m 


-  .080    -  .095    -  .115    -  .055    -  .055 

Reference  C.G.  at  0.6 
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-  .075 


2.  Pitch  Control  Law 

The  pitch  control  laws  for  the  circular  and 

elliptical  airframes  are  shown  in  Figure  2.6  as  given  in 

the  [Ref.  2].   Lag-Leads  were  used  to  prevent  guidance 

noise  saturation  problems.   An  integrator  was  used  in  the 

acceleration  error  branch  of  the  circular  control  law  to 

satisfy  the  guidance  requirement  of  zero  steady  state 

error.   In  the  elliptical  control  law,  an  integrator  was 

placed  in  the  rate  error  path  (i.e  path  after  the  adder 

(q+y) /  as  shown  in  Figure  2.6)  for  reduction  of  the 

acceleration  response  overshoot  below  10  percent  and  a  gain 

K=1.1053  was  placed  in  series  with  the  acceleration  command 

for  zero  steady  state  acceleration  error. 

A  normal  acceleration  command  n_  equal  to  1  gee, 

c 
is  applied  to  the  pitch  control  law  which  uses  measurements 

of  missile  body  pitch  angular  rate  (q)  and  pitch  normal 

acceleration  (n_)  to  determine  the  required  actuator 

command  (5   ).  The  actuator  is  modeled  as  a  first  order  lag 

pc 
at  188.4  rad/sec  and  is  shown  in  Figure  2.5. 

3.  Transfer  Functions  of  Aerodynamic  Models  for 
Circular  Airframe 

The  equations  which  are  represented  by  the  block 
diagram  of  the  aerodynamic  model,  of  the  uncoupled  Pitch 
channel  (Figure  2.5)  are: 
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F           J                4*  (II. C. 3-1) 

^^-^[srtC^Jp+C^^j/y/      [gccj]  (ii. c.  3-2) 

OL=CT.3[C|  +  3a.X^/v]     [jcf]  (II.C3-3) 

Utilizing    Laplace   Transformation,    the  above 


equation   becomes: 

°"  T 

lyy 

Hi=-^  (5?.3C^P  +  C^.a)/W 

5a  =  5T.3(c|  +  32.ai1i/V) 


(II. C. 3-4) 

(II. C. 3-2) 
(II. C. 3-5) 


Introducing    equation    II. C. 3-2    to    II. C. 3-5    the    last 
equation   yields. 


S<|  = 


S  a  -  5*.  5 


USA  (C*.<x  +  5*3  W  -fp  ) 

i * 1 i_      (II. C. 3- 

IYy 

311   /       9  &  (51-3  W fy  ^^Hq'o)  \ 


6) 


1+ 


w 


(II. C. 3-7) 
Rearranging    above   equations,    in   terms   of   q   and   a, 

we   obtain: 

5sU     -    \  /win    e*j      _       \    r 

(II. C. 3-8) 


-5"^(^^^v=(-^^wPyp(II.c.3-9, 


w     ;      v      w        wp 

Applying  the  Cramer's  rule  in  the  system  of 
equation  II. C. 3-8  and  II. C. 3-9  the  transfer  functions  q/ 6. 

and  a/6   are  obtained. 
P 
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x2 


Ivy  y         \  VMVY 

w   y         iYy 


(II. C. 3-10) 


And: 


a 


Where: 


^•S^Wp    |5   -J   3233.29  q-  ^-ci  -CmX, 

W        /        \  ivy 


P  ,2    ,     /  K-q^W   \    ^  5T-3  9  ^d  C*1a 


[Jii^s 


Lvy 


(II. C. 3-11) 


=    (57.3) (32.2)       =    0<48 


V  2 

S  =    3.14159    (ftz) 

d  =2  (ft) 

q  =  1650  (lbs/ft2) 

Iyy  =  804  (slug-ft2) 

W  =  2525  (lbs) 


C     =  -0.06 
m 
a 


Cm    =  -0.08 
m 

CN  p   =  0.15 
a 


CM     =0.04 
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a     =  0 
e 

As  referred  to  in  Appendix  A  and  Table  I.   After 
substituting  the  aforementioned  values  into  equation 
II. C. 3-10/  it  becomes: 

Introducing  equation  II. C.  3-12  in'co  equation 
II. C. 3-2  the  last  equation  yields: 


\      w      /  ~         jyy 

(II.  C. 3-13) 
Substituting  the  above  mentioned  values  into 
equation  II. C. 3-13  and  minor  manipulating,  it  yields: 

*L  r        106-0?^  x<oV- 19. sa      fwts/mrf) 

dP      22.  5^5 XiO'5-  5  V 3.3633X40 *S  +  i         (II. C.  3-14) 
The  equations  II. C. 3-12  and  II. C. 3-14  are  the 

pitch  aerodynamic  transfer  functions  for  the  pitch  angular 

rate  q  about  the  YB  and  the  achieved  maneuver  acceleration 

in  the  zD  direction  (n„) . 

a  Z 

4.   Design  Approach  and  Analysis  For  Circular  Airframe 
A  Fixed  Flight  condition  was  selected  for  these 
preliminary  performance  studies  of  circular  and  elliptical 
airframes.   The  selected  flight  condition,  at  60KFT 
altitude  and  mach  3.95,  provides  a  sufficiently  low  dynamic 
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pressure,  so  that  missile  maneuvers  will  result  in  large 
enough  angles-of-attack,  to  exercise  side  slip  control. 
The  above  fixed  condition  was  selected,  in  order  to  reduce 
the  complexities  of  design  of  the  CBTT  autopilots. 

The  uncoupled  autopilot  design  technique  was 
classical,  using  a  combination  of  frequency  response  and 
root  locus  technique  [Ref.  2]  to  achieve  practical 
bandwidths  (i.e.,  sufficient  high  frequency  attenuation) 
and  therefore  provides  the  range  of  required  missile  body 
angular  rates  and  control  motions.   The  analysis  of  the 
CBTT  autopilot  channels  is  based,  on  the  time  responses  of 
manuever  plane  accelerartion,  body  angular  rates  and  tail 
incidence  angles. 

For  purposes  of  analysis,  a  CSMP  program  was 
written  (Appendix  D)  using  the  equations  which  are 
represented  by  the  block  diagrams  of  aerodynamic  model  and 
linear  uncoupled  pitch  control  law,  as  are  shown  in  figures 
2.4  and  2.6. 

The  equations  of  the  uncoupled  linear  pitch 
aerodynamic  model  and  control  law  are: 

a.   Equations  of  Control  Law 

(1)   Acceleration  Filter  Equation 

~£o"  +  A  (II. C. 4-1) 

Rearranging  above  equation  II. C. 4-1  one  obtains: 
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X  z  -150X+  450  ^  (i i.e. 4-2) 

(2)  Acceleration  Compensator  Equation 

V= —[Vli-X)       (II. C. 4-3) 

5-(4-+i) 

Utilizing  inverse  Laplace  Transformation  and 
rearranging  the  above  equation,  it  becomes: 

V  +  5Y  =  -i.355Uii2c-X)-193.5X^63-(Tlic-X>)        (ii.c.4-4) 

Using    state-space    representation   of   a    system, 
in   which   the    forcing    function    involves   derivative   terms 
[Ref.    5:    pp.    675-678]    one   obtains: 

Y<   =  Y2  -1.3534  h^-x)  (ii. c. 4-5) 

Y2  r-5Y2 -l-G.5^2 -(Tl^-X)  (ii.c.4-6) 

Y    -    ^l  (II. C. 4-7) 

(3)  Rate  Compensator  Equation 

**"",(7**1)   m) (ii-c-4-8) 

Using  Inverse  Laplace  Transformation,  it 

y  i elds :  •  st    • 

&P,Z  -O.M3-Spc-3H0.5SX40  (Yt-oO-S.^O&'CY^XO 

c  u  (II.C.4-8a) 

b.        Actuator    Equation    (Figure    2.5): 

cjp-  - .<5pc  (II. C. 4-9) 

-hi 
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or 

£p=   -12S4  <Jp+  183.4  Sfc  (II. C. 4-10) 

c.      Equations    of   Aerodynamic   Model 

In   order    to    achieve    the    plots    of    the    time 

response    of    q,    n    ,    6      in    the    linear    uncoupled    pitch 

z    p 

channel,  the  transfer  functions  which  are  represented  by 
the  aerodynamic  model  and  derived  in  section  II. C. 3,  were 
used : 

q      _      -1.J3G1-S    -0.159  (ii. c. 4-11) 


<*P  25.5i|5Xl6V  +  3.3635X^63S  +  1 


or 

-3  o     -5 


<}  hS5HSXWSZ+  3.3633  XlO^SH-i)  z  (-1.13CJ  -5  -  0.159>-  Sp 

(II.C.4-lla) 

Rearranging  and  minor  manipulating  equation 
II.C.4-lla  it  yields: 

*<q  +  ^g.AS2Xl63q  +  Mq.3556  ^  ='59.2635  £p  -J.052-Jp     (ii.c.4-12) 
Utilizing  state  representation  of  a  system  in 

which  the  forcing  function  involves  derivative  terms,  as  in 

section  II. C. 4,  one  obtains: 

X<  -  X2  "  59.2*35- Sp  (II. C. 4-13) 

)(2    Z   -^.3556X^^49. mx*03*2+J.m£p  (II. C 4-14) 

q  =  Xi  (ii. c. 4-15) 

Also,    rearranging    and  manipulating  in    the    same 
way,    equation    II. C. 3-12   yields: 

Z|     -    12  --?OJ.QOi.\l63.  <Jp  (II. C. 4-16) 

i2   Z-m.555(,li-  149.^2X46%  -1<W  33X40   Sf>  (II. C. 4-17) 
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\-~    (*1  +  M.T05  Sp)  (ii.  c.  4-18) 

Utilizing  state-space  representation,  all  the 
aforementioned  equations  can  be  modeled  in  the  following 
ninth  order  system: 


X 


Vi 

4 
4 


o        1        o        o    o        o 

-1W/5M  -0.15      O  O    O  O 

0        0        O        i    o       O 


0 

O 

0 
o 


0     -*W.34   -0,45  0        O 
O    +150       0    150      O 


O 
O 


X23i    OMi 
O         O 


0 
O 
O 
O 


O 
O 
O 

o 
o 

i 


O  -59.36 

O  1.T8I5 

O  -0.302 

0  -xs4s 


O      <M05 

4 


0    O     1.353 
0  0     6.5?a    -5 


[V 

0 

Xi 

0 

?i 

0 

1 

?2 

0 

•  X 

4 

0  i 

Yi 

-135 

Y2 

CSf 

4 

a4u 

Jp 

0 

o       o 

O        O 
0    0     -0.46f  -0.3W1  -0.^3  -30.2 

o  o      o       o     l%%4  -i$M_ 

(II. C. 4-19) 
where  the  state  variables  are: 

pitch  angular  rate  state  variables 
(section  II. C. 4) 

Achieved  Body  Acceleration  state 
variables  (section  II. C. 4) 
X      :   output  of  Acceleration  filter 

6P 

c    :   input  command  in  actuator  network 

Sp     :   tail  incidence  angle 

The  state  X,  X,,  Y,  ,  6   of  control  law  are  shown  in 

c 
figure  2.6. 


Xl'  X2 


Yl'  Y2 
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Executing  the  program  of  Appendix  D,  using  an 
input  step  function  which  represents  the  "lGee  Command"  and 
trim  angle  of  attack  a  =0,  the  time  response  plots  of 
pitch  normal  acceleration  nz,  pitch  angular  rate  q  and 
pitch  tail  incidence  6   for  circular  airframe  are  obtained. 

Figure  2.7  shows  the  pitch  acceleration  response 
of  the  circular  airframe  due  to  a  one  Gee  acceleration 
command,  when  the  Missile  Aerodynamic  Models  are  linearized 
about  zero  angle  of  attack  (a) .   The  response  has  a  0.5 
seconds  time  constant,  steady  state  error  0.005  and  2.5% 
overshoot,  which  are  according  to  the  requirements  referred 
to  in  appendix  C. 

Figures  2.8  and  2.9  show  the  required  body 
angular  rate  and  control  surface  deflection,  to  achieve  the 
acceleration  response.   These  responses  matched  those 
presented  in  [Ref  2]. 

5.   Design  Approach  and  Analysis  For  Elliptical 

Airframe 

The  design  technique  used  for  the  elliptical 
airframe,  is  the  same  as  that  of  the  circular  airframe. 
The  analysis  of  the  Pitch  channel  for  the  elliptical 
airframe  was  based  again,  on  the  time  response  plots. 

For  the  above  purposes,  a  CSMP  program  was 

written,  as  in  the  case  of  circular  airframe,  referred  to 

in  Appendix  D. 
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Fig.  2.7   Pitch  Normal  Acceleration  (nz)  vs  Time  (t) 
Uncoupled  Pitch  Channel 
Circular  airframe  (lGee  Command,  a   =0). 
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2. 


n 
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Fig.   2.9   Pitch  Tail  Incidence  (5  )  vs  Time  (t) 

Uncoupled  Pitch  Channel 

Circular  airframe  (lGee  Command,  ae  =0) . 
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The  equations  of  the  uncoupled  linear  Pitch 
channel  for  elliptical  airframe,  which  are  used  in  the 
above  CSMP  program,  are: 

a.   Equation  of  Control  Law 

The  control  law  of  the  uncoupled  linear  Pitch 
channel  for  the  elliptical  airframe  is  very  similar  to  the 
control  law  for  the  circular  airframe,  with  only 
differences  in  the  gain  of  acceleration  and  rate 
compensators.   Elliptical  control  law  is  shown  in  figure 
2.6.,  taken  from  [Ref.  2], 

(1)  Acceleration  Filter  Equation.   The 
Acceleration  filter  equation  is  the  same  as  in  section 
II. C. 4 

(2 )  Acceleration  Compensator  Equation 

Y=  (-O.0Os/|-  +  iy  (i.i053»lvx)  ("-C.5-D 

Utilizing,  inverse  Laplace  Transformation 

and  rearranging  II. C. 5-1  it  becomes: 

Vr-€Y-  530.5^X10^  40.M2-X    (II. C. 5-2) 


(3)   Compensator  Equation 


Y-X<  5" 

(II. C. 5-3) 
Using  inverse  Laplace  Transformation,  it 


yields : 

-3 


5pc  -  -3*3.?5X<03(y-X0  -3.0Y(Y-Xf)  (i i.e. 5-4) 
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b.  Actuator  Equation  (Figure  2.5) . 
As  in  section  II.C.4fb 

c.  Equation  of  Uncoupled  Linear  Pitch 
Aerodynamic  Model 

The  equation  for  pitch  angular  rate  was 
derived,  as  in  section  II. C. 3,  the  only  differences  are  in 
the  constants  W=2475(lbs), 
lyy=790 (Slug-ft2) ,  Cm  =0.015,  Cm  =-0.055,  Cn  =0.02  and 

a  6  <S 

P  P 

C   =0.18,  as  referred  in  the  Appendix  A  and  Table  I. 
a 

Manipulating  according  to  the  inverse  Laplace 

Transformation  rules  and  rearranging,  one  abtains. 

~-=(3.66S+0.6SS)/(-88.^Xi6V+^5.95X^d3S-H)       (II. C. 5-5) 

Using  state-space  representation  of  a  system, 
in  which  forcing  function  involves  derivative  terms  [Ref. 
5:  pp.  675-678.]   It  yields: 

Xi  =:  X2  -4i.29-£p  (II. C 5-7) 

X2  z  ii.2«5-X4+4Tg.Q4x465x2-i5.*g^-<5'p  (ii.c.5-8) 

S  =  *1  (II.C.5-8a) 

Also,    Manipulating    with    same    way    the    equation: 

^      ,    -212.4jiXl6352+gi.6  (ii. c.  5-9) 

Jp  -8UMX^S2+15-95XH0*S+  A 


Yielding 


Z,    =   22  +  52*.3X10.<SV  (II. c. 5-10) 

i2   =  11.2815-^  +  47944  X40^2- W* 35  ^P  (II. C. 5-11) 

Y\2  r   -  (i-i  +  2-936 -5f)  (II. C 5-12) 
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Utilizing  state-space  representation,  all  of 
the  above  equation  can  be  modeled  in  the  following 
eighth-order  system: 


5P 


l 


o 


o  o   -ii-« 


1J.28  CIS 

0 

0 

0 

o    o     -15.ft 

O         0 

o 

i 

0 

O  O     0.5233 

0          0 

iJ.28 

OAZ 

0 

0    0    -8W.3 

0         0 

150 

0 

H50 

0    0     »W05 

O         0 

0 

0 

0.M3 

-4  0       o 

3.o*  aw 

o 

o  - 

-0.W 

0U8  0     HS.W 

0          0 

0 

0 

0 

0   #14  -<«4| 

Xi 

0 

Xz 

0 

2) 

0 

5' 

0 

Y 

•0.53 

ilk 

-0.203 

Sp 

0 

^r 


(II. C. 5-13) 
where  the  state  variables  are: 

X,,  x~:   pitch  angular  rate  state  variables 

(section  II. C. 4) 
Y,,  Y2 :   Achieved  Body  Acceleration  state 

variables  (section  II. C. 4) 
X      :   output  of  Acceleration  filter 

c    :   input  command  in  actuator  network 

6p     :   tail  incidence  angle 

The  state  X,  X-,  ,  Y,,  5  of  control  law  are  shown  in 

c 
Figure  2.6 
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Utilizing,  the  above  state  variables  system,  a 
CSMP  program  was  written  (Appendix  D)  to  achieve  the  time 
response  plots  of  n  ,  q  and  5  .   As  input  a  step  function 
is  used,  which  represents  a  one  gee  command.   The  trim 
angle  of  attack  (a  )  is  equal  to  zero. 

Figure  2.10  shows  the  pitch  acceleration 
response  of  the  elliptical  airframe  due  to  the  command.  The 
response  has  a  time  constant  of  about  0.5,  steady  state 
equal  to  zero,  and  overshoot  o% ,  keeping  the  requirements 
referred  to  in  Appendix  C. 

Figures  2.  Hand  2.12  show  the  required  body 
angular  rate  and  control  surface  deflection,  to  achieve  the 
acceleration  response. 

6.   Comparison  of  Airframe  Response 

Figures  2.7  and  2.10  show  that  the  Pitch 
acceleration  response  of  the  circular  and  elliptical 
airframe  are  approximately  the  same  having  a  0.5  second 
time  constant  and  negligible  overshoots. 

Figures  2.8,2.9,2.11  and  2.12  show  that  to  achieve 
the  acceleration  response  the  elliptical  airframe  requires 
less  body  angular  rate  and  control  surface  deflection.   The 
elliptical  airframe  requires  less  control  surface 
incidences,  because  the  transfer  functions  have  a  higher  DC 
gain  and  poles  which  are  located  at  a  lower  frequency  due 
mainly  to  a  more  neutrally  stable  airframe. 
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Fig.  2.10   Pitch  Normal  Acceleration  (nz)  vs  Time  (t) 
Uncoupled  Pitch  Channel 
Elliptical  airframe  (lGee  Command,  ae=0) 
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Fig.  2.11   Pitch  angular  rate  (q)  vs  Time  (t) 
Uncoupled  Pitch  Channel 
Elliptical  airframe  (lGee  Command,  a  =0) 
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Fig.  2.12   Pitch  tail  incidence  (5  )  vs  Time  (t) 
Uncoupled  Pitch  ChannelP 
Elliptical  airframe  (lGee  Command,  a  =0) 
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The  DC  gain  [(qS/w).(CN   .Cm   /   Cm  -  C^    ) ]  is  directly 

a      6         a      fi 
P  P 

proportional  to  pitch  control  moment  C  and  inversely 

P 
proportional  to  magnitude  of  stability  margin  in  Pitch 

(i.e.,  C  /CM  or   distance  from  center  of  pressure  to 
m   in 
a   a 

center  of  gravity) .   Thus  the  nearly  neutrally  stable 

elliptical  airframe  is  expected  to  have  a  higher  gain  than 

the  stable  circular  airframe  even  though  its  control  moment 

coefficient  is  slightly  smaller.   Higher  DC  gain  will 

require  less  control  surface  incidence.   The  zeros  of  n„/$ 

z  y 

are  directly  proportional  to  C    and  the  location  of  the 

a 
ratio  C    /  C     or  the  distance  from  the  point  of 
6        6 
P       P 

action  of  tail  forces  to  the  center  of  gravity. 

Hence,  for  the  elliptical  airframe  which  has 

larger  CN   the  zeros  of  n  /$   are  located  at  a  higher 

frequency.   Lower  control  surface  incidences  and  higher 

frequency  n  /<S   zeros  of  the  elliptical  airframe  will 
z   p 

simplify  the  design  of  a  rapidly  responding  Pitch 
autopilot. 


D.    UNCOUPLED  LINEAR  YAW  CHANNEL 

The  purpose  of  the  Yaw  channel  autopilot  of  a  CBTT 
autopilot  is  to  minimize  sideslip  angle  or  provide 
coordinated  motion  between  Roll  and  Yaw  channels.   This  was 
accomplished  in  two  ways.   First,  the  uncoupled  Yaw  channel 
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autopilot,  (i.e,  Roll  and  Pitch  dynamic  effect  neglected) 

was  designed  as  a  regulator  (i.e  no  guidance  command  and 

with  rate  and  acceleration  feedback)  to  help  minimize 

sideslip  angle.   Second,  to  aid  in  sideslip  control,  the 

regulator  was  commanded  from  the  Roll  channel  as  explained 

in  chapter  III. 

1.   Aerodynamic  Model  Control  Law 

A  block  diagram  of  the  uncoupled  Yaw  channel  is 

shown  in  Figure  2.13,  taken  from  [Ref.  2].   In  this  diagram 

the  aerodynamic  model  and  the  Yaw  control  law  are  involved. 

Also,  the  Yaw  control  law  without  the  accelerator  is  shown 

in  Figure  2.14.   These  figures  were  abstracted  from 

[Ref.  2]. 

The   normal  acceleration  n     is  not  used  to 

yc 
command  the  CBTT  autopilot.   Instead,  it  is  used  for  the 

design  and  analysis  of  the  uncoupled  channel.   The  command 

used  by  the  coupled  system  is  shown  in  dashed  lines  and  is 

a  yaw  angular  rate  command,  r  . 

The  yaw  control  law  is  governed  by  missile  body 

yaw  angular  rate  (r)  and  yaw  normal  acceleration  n  .   The 

yaw  control  law  determines  the  required  command  6 

yc 

to  an  actuator,  which  is  approximated  as  a  first  order  lag 
at  30  Hz  (188.4  rad/sec) . 

The  aerodynamics,  linearized  about  a  trim  a 
angle-of-attack  equal  to  zero. 
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2.   Transfer  Functions  of  Aerodynamic  Model 

The  aerodynamic  characteristic  of  the  Yaw  channel 
are  modeled  by  a  second-order  system,  by  the  following 
equation,  where  the  states  are  r,  yaw  angle  and  3  , 
sideslip  angle: 

T=ft£lflRWGv«  +  9ttCaSrSv)   Cvocol/s«c)  (ii.d.2-1) 

«:5T.3(-r+-^.*v)       Utj)  (ii.d.2-2) 

and    also,    by    the    equation: 

V\yr  (q^/W)-   (Cy  .0  +  57.3  Cyjy.Sv  )    Ipts)  (II. D. 2-3) 

substituting    equations    II. D. 2-3    into    II.D.2-2,    the 
last   equation   becomes: 

6=  5T.3[-t+^(-|f  C«r««  +  SMCrV^y)")J        (ii.d.2-4) 

Taking  the  Laplace  Transform  of  the  equations,  one 
obtains^ 

V5=  -9^-  •  (Cn6-«+W3C^  Jy)  (ii.d.2-5) 

for  S*3  (-V-+^?(^  (fy-g+ST-B   Cvfi.Sy))) 

(II. D. 2-6) 
Rearranging    above   equations,    in   terms   of   q   and  a, 
they   become: 

n--^*— <*•*  =  '"^      HJy        ui.D.2-7, 

(II. D. 2-8) 


Applying  the  Cramer's  rule  in  the  system  of  the 
equations  II. D. 2-7  and  II. D. 2-8  the  transfer  functions  r/6, 
and  B/<S     are  obtained. 

y  ,,v.9 .         v        __  ../^^\2    . 


r 

-  - 

^        i«             °V  -  '        v.  i« 

\          ~M< 

57    -i<j  ■  -"a-vrdyj 

Jy 

^(-^■^H  +  C*^-} 

cne 

(II. D. 2-9) 

and 

= 

h 

,2.    /  KqSGia  >|  ,       5T.3 q S'd Cmot 
S  +   [      W         ;  J               lyv 

(II. D. 2-10) 

K 

(57.3) . (32.3) 
—         —         _____              —    a    ah 

—                        _____                                        _          |0   B     ^Q 

V 

q 

=      1650    (lbs/ft2) 

s 

=      3.14159     (ft2 

d 

=      2    (ft) 

Izz   =   810 (slug-ft2) 

w    =   2525  (lbs) 
C„    =   -0.025 

n3 
C 

n6   =   -0.050 

C 
y6   =   -0.065 


C 
Y«   =   0.021 

a     =   0° 
e 


As  referred  to  in  Appendix  A  and  Table  I. 
Substituting  the  aforementioned  values  into  equation 
(II. D. 2-9),  it  becomes: 


~5T -5...5W53- i»x«f*s  ti (""c)  <"•»•*■«> 

Introducing    equation    II. D. 2-10    into    II. D. 2-3   and 
rearranging    the    last    one,    the    transfer    function   of    yaw 
normal    acceleration    in    terms    of    the    yaw    control    incidence, 
is    obtained: 

( '5T3Kfq^)2  r       r     1  C       33  83.39  fitfd  ,        , 


*       s'+pjL^Vs-    »»**'€ 


■Woe 


(II. D. 2-12) 

After    substituting    the    values    into    the    equation 
(II. D. 2-12),    one    obtains: 


-viy       -13.  w  s2+  IT.?6 


^TT 


(gas/moO 


Jy         -5*J.5HX10S-Z8X<0-S    +1 

(II. D. 2-13) 
The   above    transfer    functions    are    used    for    finding 
the    time    responses    of    the    yaw   angular    rate    r    and    the    yaw 
normal    acceleration. 

Using    inverse    Laplace    transformation,    the    equation 
II. D. 2-11   and    II. D. 2-13    become: 
?  4-51338X10 •f-J$.3335-r  r-36.?6l-Sy-2.i5-Sy  ( II. D.  2-14) 

My +51338Xi6i'«\ -A83S35  ">1y  =  2.38*  iy -325. 629- £v  (II.  D.  2-15) 

The    equations    II. D. 2-14    and    II. D. 2-15    involve 

derivative    terms    in    the    forcing    function.      Using    rules   of 

state-space    representation   of    a    system    [Ref.    5:    pp 

675-678] ,    they    turn: 
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Xj  =  X2-36.?6J-cW  (ii. d.  2-16) 

X2   r  18.335-Xi-51.338X'l6.3C2-g63.omx^0^y       (H.D.2-17) 


V     -    Xi  (II. D. 2-18) 

2l=    22-m5«WX463.*/  (II.D.2-19) 

?2  =  iS.335-2<i -51.338X403*2- 29J.S51oy      (II. D.  2-20) 
tyz    2i  +  2.3?f.Jy  (ii. d. 2-21) 

3.      Equations   of  Yaw  Control    Law 

As   mentioned   before,    the    yaw  control    law   is 

governed   by  missile   body   yaw   angular    rate    ( r)    and    yaw 

normal    acceleration   n      and    it    is    used    to   determine    the 

required    command    6Y      to    the    actuator.      The    yaw   control    law 

c 
at    the    flight    condition   of    interest    (i.e,    Mach    3.95,    60kft) 

is    shown    in   Figure    2.14,    taken    from    [Ref.    2].      The    rate 

compensation   determines   the    high    frequency   attenuation   and 

is    used    to  minimize    aerodynamic   variations   on    the   quality 

of    regulation.      The    acceleration   compensation   determines 

the   acceleration   bandwidth  via    the    time   constant   of    the 

acceleration   response   due    to    a   step   command   of   acceleration 


at  nY    . 


a.   Acceleration  Compensator  Equation: 

(II. D. 3-1) 


0.25+1     {    Yc       YJ 


Inverse  transforming  and  rearranging  the 
II. D. 3-1,  it  turns  into: 

V  =  -5Y +X.5913  l^y-^Yc)    (ii.d.3-2) 
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b.   Rate  Compensator  Equation: 


fr  _    H«5  C-j-5+0    (y+r)       (II. D. 3-3) 


Utilizing  inverse  Laplace  Transformation 
and  introducing  equation  II. D. 2-16,  the  equation  II. D. 3-3 
turns  into: 

<Jyc  Z   0.MS5  (Y  +  XO  +4.  85(Y  +  XO  (II.D.3-3a) 

c.   Actuator  Equation: 

Jyr  _i Jyc       (h.d.3-4) 

+  i 


183.4 

Inverse  transforming  and  rearranging  the 

equation  II. D. 3-4,  it  turns  into: 

£y  r-iSS.4-Sv  +  18S-4  lie  (ii. d. 3-5) 

Using  state-space  representation,  all  the 
aerodynamic  model  and  control  laws  can  be  modeled  in  the 
following  seventh  order  system: 


Xt 

k 

Y 
& 

Si 


o       1       o        o       o 

i*.335  -asm     O         0         o 

0        o       o        1        o 


o    -am* 
o    -m.% 

0       3.812? 


o  o    is.  335  -asm    o 

0  0      J.59*       0        "5 

H.85  0.MS5    0.7*5       0      2.425       0     -if.SM 

0  0         0          0          0     18*4  -l««-4 


Xi 

0 

x2 

o 

2i 

0 

• 

22 

+ 

0 

Y 

-i.59? 

ivc 

0-W 

*lvc 


(II. D. 3-6) 
where  the  state  variables  are: 

X,,  X2:   pitch  angular  rate  state  variables 

(section  II. D. 2) 
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Y,,  Y2:   yaw  acceleration  state  variables 

(section  II. D. 2) 
Y      :   output  of  acceleration  compensator 

c    :   input  command  in  the  actuator 

Sp     :   yaw  tail  incidence 

The  state  X, /  Y,  6y  of  control  law  are  shown  in 

c 
Figure  2.14. 

4.   Design  Approach  and  Analysis  of  Circular  Airframe 

The  design  technique  used  for  the  Yaw  channel  is 
the  same  as  the  technique  used  in  the  Pitch  channel.   The 
analysis  of  it,  for  circular  airframe  is  based,  on  the  time 
responses  of  the  yaw  angular  rate  r,  the  yaw  normal 
acceleration  n   and  yaw  tail  incidence  5  . 

For  purposes  of  analysis,  a  CSMP  program  was 
written  (Appendix  E)  using  the  state-variable  system  which 
is  represented  by  the  block  diagram  of  the  yaw  aerodynamic 
model  and  control  law. 

Executing  the  referred  program,  using  as  input  a 
step  function,  which  represents  the  "one  Gee  Command"  and 
trim  angle  of  attack  a  =0 ,  the  time  response  plots  are 
obtained . 

Figure  2.15  shows  the  yaw  acceleration  response  of 
the  circular  airframe.   The  response  has  a  0.39  seconds 
time  constant,  7%  overshoot  and  steady  state  error  0.018 
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Fig.  2.15   Yaw  normal  acceleration  (n  )  vs  time  (t) 
Uncoupled  Yaw  Channel 
Elliptical  airframe  (lGee  Command,  ae=0) 
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which  are  according  to  the  requirements  referred  in 
Appendix  C. 

Figures  2.16  and  2.17  show  the  required  body  yaw 
angular  rate  and  control  surface  deflection,  to  achieve  the 
acceleration  response.   Again  these  results  matched  those 
of  [Ref.  2]  . 

5.   Aerodynamic  Transfer  Function  for   Elliptical 

Airframe 

The    transfer    functions    for    yaw   angular    rate    and 
yaw   normal    acceleration    for    the    elliptical    airframe,    were 
derived    in    the    same   way    as    in    section    II. D. 2.      The 
only   differences    are    in    the    following    constants 

w  =    2475     (lbs) 

Izz        =    853    (slug-ft2 

C  =    0.024 

C  =    -0.042 

6 
Y 

C  =    -0.043 

B 

CY  =    0.016 

6 
Y 

As    referred    to    in   Appendix   A  and    table    I. 
a.      Transfer    function   of    yaw   angular    rate   r. 


r     .      -Ltt.S-S-0.OS9l  (^(ii.d.,-1, 

6^  59.838X10    SV  2.593X10 -S  +  i 
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Uncoupled  Yaw  Channel 
Circular  airframe  (lGee  Command,  ae=0) 


67 


8 

a 

& 

o 

8 

o 

S 

a 

s 

o 

o 

o 

§8 

5o 

a  _ 

a 

o 

8 

o 

5 

i 

o 
1 

8 

i 

? 

8 

? 

1/ 

S 

v 

0.00  0.2  0.50  0.75   1.00   2 .25   1.S0   1.79  2.00  2.25  2.50  2.75  3.00  3.25 

TCSEC1 

Fig.  2.17   Yaw  tail  incidence  (6y)  vs  Time  (t) 

Uncoupled  Yaw  Channel 

Circular  airframe  (lGee  Command,  ae=0 ) 


68 


b.   Transfer  function  of  yaw  normal  acceleration  nY 

Ik--  mMixto'?-  U (*SL\iu  d  5-2 

Si  59.  S3S  X  i  0~3.  S2+ 2.593X4 6*  S+i    V  tad  ) 


Using    inverse    Laplace    transformation,    the   equation 
II. D. 5-1    and    II.D.5-2    becomes: 

V+U333XW?f4*mM  =-29.2«.£y- 9^.36X^5?  JV       ("-D.5-3) 

^Y-f43.333X^63'yiy4^?Ji?-7ly==^qi5?-<Jy-^8-653-(5,y       (II. D. 5-4) 
Utilizing    state    representation   of   a    system    in 
which   the    forcing    functions    involve   derivative    terms,    as    in 
section    II. D. 2,    the    above    equations    yield: 

X<   -    X2-  29.221  •  <JV  (II. D. 5-5) 

Sfe  =  -iUmX^-M3.3^Xl6^X2+266.OTX^(J,y      (ii.d.5-6) 

f:    Xl  (II. D. 5-7) 

4  -    22-  93.013X40. <5v  (ii. d. 5-3) 

22"=  H6.W:J.£i-43.333Xl6?22-'150.W3?-5y  (ii.d.5-9) 

^Y~   2*4-  i.95i*-0Y  (II. D. 5-10) 

The  aerodynamic  model  for  elliptical  airframe  is 

shown  in  Figure  2.13. 

6.   Yaw  Control  Equation  for  Elliptical  Airframe 

The  control  law  of  uncoupled  linear  Yaw  channel  for 
elliptical  airframe,  is  similar  to  the  control  law  for 
circular  airframe.   There  are  differences  in  the  gains  of 
the  acceleration  and  rate  compensators,  as  is  shown  in 
Figure  2.14. 
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a.  Acceleration  Compensator  Equation 

V=  (0.S393/O.25-S+0*  (-'«YC+'V)         (II. D. 6-1) 

Utilizing  inverse  laplace  transformation  and 
rearranging  the  equation  II. D. 6-1,  it  turns  into: 

Y  =  -4V  +  3.35?4  ("nyc  +  ^y)         (ii.d.6-2) 

b.  Rate  Compensator  Equation 

SYt.  6.08- Ofc+0    (Y_r, 

S  (II. D. 6-3) 

Inverse  Laplace  transforming  and  introducing 
the  equation  (II. D. 5-7)  into  the  equation  (II. D. 6-3),  it 
turns  into: 

<$YC  =  0-608  (Y+X<)  +  6.0$(Y+XA       (ii.d.6-4) 

c.  Actuator  Equation 

Same  as  the  equation  II. D. 3-5. 
Using  state-space  representation  all  equations  of  the 
aerodynamic  model  and  control  law,  can  be  modeled  in  the 
follow  seventh-order  system: 


x, 
i, 

» 

4 

5y 


0         1 

O 

0 

o 

-iG.32  -0.U33 

0 

0 

0 

O          O 

0 

1 

0 

0         O      -i6.?2  -0.435      O 
0         0       3.56       0       "I 
6.08   0.608  5.041       0     '3M5Z 


O  -29.22 

O  0.26? 

o  -cm 

0  -<50.6 

0  6.553 

0  3.<&<| 


O      1*8.<f    'iW 


x2 

^2 
Y 


\ 


o 
o 

0 

0 

-3.36 

-104 

0 

(II. D. 6-5) 
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where  the  state  variables  are: 

X, ,  X^ :   yaw  angular  rate  state  variables 

(section  II. D. 5) 

Z, ,  Z«:   yaw  normal  Acceleration  state 

variables  (section  II. D. 5) 

Y      :   output  of  acceleration  compensator  network 

6y    :   input  command  in  the  actuator 

c 
6      :   yaw  tail  incidence 

The  state  X-,,  Y,  <5y  °f  control  law  are  shown  in  Figure 

c 
2.14. 

7.   Design  Approach  and  Analysis  for  Elliptical 

Airframe 

The  design  technique  and  analysis  for  elliptical 
airframe  is  the  same  as  that,  in  section  II. D. 4.   Again, 
for  purposes  of  analysis,  a  CSMP  program  was  written 
(Appendix  E)  using  equations  of  aerodynamic  model  and 
control  law. 

Figure  2.18  shows  the  yaw  acceleration  response  of 
the  elliptical  airframe.   The  response  has  0.35  second  time 
constant,  overshoot  about  6%  and  steady  state  error  0.01 
which  are  according  to  the  requirerments ,  (Appendix  C) . 

Figures  2.19  and  2.20  show  the  required  body  yaw 
angular  rate  and  control  surface  deflection,  to  achieve  the 
acceleration  response. 
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E.    LINEAR  UNCOUPLED  ROLL  CHANNEL 

Generally,  the  purpose  of  a  roll  channel  autopilot  is 
to  command  the  <t>  signal,  in  order  to  roll  the  missile 
through  an  angle  <t>,  measured  from  the  vertical  axis. 

The  roll  channel  of  a  coordinated  bank-to-turn 
autopilot  is  commanded  to  roll  the  missile  so  as  to  put  the 
preferred  maneuver  direction  of  the  missile  in  the 
direction  of  the  guidance  acceleration  command,  while  the 
pitch  channel  acceleration  is  commanded  to  produce  the 
total  magnitude  of  the  guidance  acceleration  command. 

The  desired  maneuver  plane  acceleration  should  be 
attained  as  rapidly  as  the  achieved  body-fixed  pitch 
acceleration.   To  accomplish  this,  the  uncoupled  roll 
channel  autopilot  was  designed  to  have  the  roll  angle  time 
constant  equal  to  the  time  constant  of  the  normal 
acceleration  achieved  by  the  uncoupled  pitch  channel 
autopilot.   This  is  not  necessarily  the  optimal 
relationship  betweem  the  two  time  constants. 

1.   Roll  Aerodynamic  Model  and  Control  Law 

A  block  diagram  of  the  uncoupled  Roll  channel  is 

shown  in  Figure  2. 21, taken  from  [Ref.  2].   In  this  diagram 

the  aerodynamic  model  is  involved.   Also,  the  roll  control 

law  is  shown  in  Figure  2.22.   These  figures  were  taken  from 

[Ref.  2]. 

The  roll  control  law  is  commanded  by  roll  angle  <t>c 

is  governed  by  roll  angular  rate  p  and  roll  angle  <j)  and 
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determines  the  required  command  5R   to  an  actuator  which  is 

c 
approximated  as  a  first  order  lag  at  30Hz  (188.4  rad/sec) . 

The  aerodynamics,  linearized  about  a  trim  angle-of-attack 

(a  )  equal  to  zero. 

2.  Equations  of  Aerodynamic  Model  for  Circular 
Airframe 

The  aerodynamic  model  is  represented  by  the 
following  equations: 

p=    (q^d/locyC^.fo  (radku1)  (ii.".2-l) 
<&=     PCradlsec)  (II-E-2"2) 

where : 

q  =  1650  (lbs/ft2) 

s  =  3.14159  (ft2) 

d  =  2  (ft) 

Ixx  =  40  (slug-ft2) 

C^  =  0.031 

As  referred  to  in  Appendix  A  and  Table  I. 
Substituting  the  above  values  in  equation  II. E. 2-1,  it 
yields: 

P  r  g.03^6  .$$  (II. E. 2-3) 

3.  Equations  of  Aerodynamic  Model  for  Elliptical 
Airframe 

The  aerodynamic  model  equations  are  similar  to 
those  for  the  circular  airframe,  with  the  only  differences 
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in    the    constant    I      =110 (slug-ft    )    and    C        =   0.023. 

XX  Jo 

fiR 

Hence    they   are: 

p  -    a.^.Jp  Crctdlsec2}  (ii.e.3-d 

<j>  =     P  OradlseO  (ii.E.3-2) 

4.   Equations  of  Control  Law  for  Circular  and 
Elliptical  Airframe 

The  yaw  control  law  for  circular  and  elliptical 
airframe  is  shown  in  Figure  2.22. 

a.   Roll  Angle  Compensator  Equation 
2. 2 


X3 


.  (cj>c-4^ 


+i  (II. E. 4-1) 


9 

inverse  Laplace  transforming  the  equation  II. E. 4-1 

it  becomes: 

X  "  -*X-H?«6  (*fc-<t>)  (II. E. 4-2) 

b.  Rate  Compensator  Equation 

v/_    0.1St»S.(JL+Q      ^  (II.E.4_3) 

Utilizing  inverse  Laplace  transformation,  and 
minor  manipulating  the  equation  II. E. 4-3,  one  obtains: 

^  r-5-Y+ 50.333X40 •(*-£)  +  0.1565 (x-p)     (ii.e.4-4) 

c.  Psuedo  Differentiator  Equation 

v  O.O1303  ^ 

Xj-      .    P  (II.  E.  4-5) 

Using    inverse    Laplace    transformation   and 
rearranging    the   equation   II. E. 4-5,    it   becomes: 

X<  r  -6Xjj  4-0.01S2-?  (ii. e. 4-6) 
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d.  Equation  of  Actuator  Compensator 

This  network  of  roll  control  law  is  different 
for  circular  and  elliptical  airframe. 

(1)  Circular  Airframe 

Using,  inverse  Laplace  transformation  and 
rearranging  the  equation  II. E. 4-7,  it  yields: 

<Vc  =  "H«>SBc+CM5<i5(H)  +  15(,'-X,>   (II. E. 4-8) 

(2)  Elliptical  Airframe 

Inverse  Laplace  transforming  and 
rearranging  the  above  equation,  it  becomes: 

« 

(II. E. 4-10) 

e.  Actuator  Equation 

5"Rc  =  ^"^ ^C  ^^«j)     (II. E. 4-11) 

Utilizing  inverse  Laplace  transformation,  the 
equation  II. E. 4-11  becomes: 

$'%z   -«M-$*  +  i0*95/52£fc(*d*ci(II.E.4_12) 

Utilizing  state-space  representation,  all  the 
aforementioned  equation  of  the  aerodynamic  models  and 
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control  laws  for  the  circular  and  elliptical  airframes,  can 
be  modelled  in  the  state  variables  systems  given  below: 


* 

x 

\ 

Sic 


X 

i 

Y 
X, 

Sit 


(1)       For  the    circular    airframe 

O         0  O         0  O         O      8.05MS 

i          0  0          0  O          O         O 

0        H.6  '2         O  0          0          0 

-0.156  0.1%  0-5  0      0.353  -<W0»W 

0        0  0         O  -6        O      6.283 

-0.103  0-U1  0.4*3   M.549  -iW.lfi    '5     -0-03? 

0         O  O         O  o      10?95.32  -ig? -^ 


F> 

O 

<D 

o 

X 

-iU 

• 

Y 

+ 

-0.8*6 

Xi 

0 

&c 

-0.121 

5, 

O 

(I 

C.E.4- 

<t>< 


o 
1 


o 
0 


(2) 
O 

O 


For    the    elliptical    airframe 


O 
0 


O 

O 


0         0 


0       tt.6     -J 
-O.ttfc  0$%      0 


OOOO 
-5      0.353     0     -0.mwn 

0-6         0     6.2S3 


-W  -ason  0.203  -3.M  59.  W3  -15  -0.58* 
0  0         0         O  O     10T95.3  -fl»*J 


p 

0 

* 

0 

X 

-iU 

•  Y 

+ 

-0.356 

Xi 

0 

S*c 

-0.504 

<J* 

o 

<*>c 


(II. E. 4-14) 


where  the  state  variables  are: 
P  :  roll  angular  rate 
<{>   :  roll  angle 

X   :  output  of  roll  angle  compensator 
Y   :  output  of  rate  compensator  network 
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X,  :  output  of  pseudo  differentiator 
6D  :  input  command  in  the  actuator 

K 
C 

5R    :      roll    tail    incidence 

The    state    4»,X,P,Y,X1,    6y      are    shown    in    Figure 

c 
2.22. 

5.   Design  Approach  and  Analysis  of  Uncoupled  Linear 

Roll  Channel  for  Circular  and  Elliptical  Airframe 

The  design  technique  and  analysis  is  the  same  as 

that  which  was  used  for  the  Pitch  and  Yaw  channel. 

For  purposes  of  analysis,  two  CSMP  program  were 

written  (Appendix  F)  using  the  state  variable  system 

represented  by  the  aerodynamic  model  and  control  law, 

Block. diagram. 

Figures  2.23  shows  the  roll  angle  response  of  both 

elliptical  and  circular  airframe.   The  time  constant  of  the 

roll  angle  response  is  0.55  seconds,  the  overshoot  3%  and 

the  steady  state  roll  angle  error  equal  to  zero. 

Figure  2.24  shows  the  required  body  roll  angular 

rate  for  both  elliptical  and  circular  airframes.   Only  the 

roll  tail  angular  deflection  figures  2.25  and  2.26  is 

different  for  the  airframes,  due  to  the  method  of 

compensating  for  a  reduction  in  the  elliptical  aerodynamic 

roll  gain.   This  is  desirable  for  the  aerodynamic  roll  gain 

(57.3qsdC4   /JXx^  to  be  as  lar<3e  as  possible  to  minimize 
6R 
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control  surface  motion.   It  is  also  desirable  to  have  as 
large  C    as  possible  to  minimize  the  effects  of 
aerodynamic  control  cross  coupling.   The  circular  airframe 
has  a  considerably  larger  aerodynamic  roll  gain  due  to  a 
much  smaller  roll  inertia  and  a  larger  control  derivative 
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III.   LINEAR  DESIGN  AND  ANALYSIS  OF  COUPLED  AUTOPILOTS  FOR 
CIRCULAR  AND  ELLIPTICAL  AIRFRAMES 

A.  GENERAL 

The  linear  design  and  analysis  technique  began  with 
uncoupled  autopilot  channels.   The  uncoupled  autopilot 
design  technique  was  classical,  using  a  combination  of 
frequency  response  and  root  locus  techniques  [Ref.  2],  to 
achieve  practical  bandwidths  (i.e,  sufficient  high 
frequency  attenuation)  and  in  turn  provides  the  range  of 
required  missile  body  angular  rates  and  control  motions,  as 
mentioned  in  section  II. 

In  this  section,  CBTT  control  laws  are  analyzed  to  add 
control  coupling  for  coordinated  missile  motion.   A  measure 
of  sideslip  control  is  obtained  by  applying  a  roll  angle 
command  to  the  linearized  CBTT  autopilot.   The  relative 
stability  of  the  autopilot  branches  and  means  for  improving 
stability  are  discussed.   An  examination  of  the  autopilot 
sensitivity  to  aerodynamic  cross-coupling  is  made. 

B.  AERODYNAMIC  MODELS  FOR  THE  CIRCULAR  AIRFRAME 

The  aerodynamic  models  of  pitch  and  roll-yaw  channels 
for  the  circular  airframe  are  shown  in  Figures  3.1  and  3.2, 
taken  from  [Ref.  2].   The  models  are  linearized  and  were 
developed  for  stability  studies  in  the  frequency  domain. 
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The  assumptions  for  linearization  are  the  same,  as  those  in 
section  II.C.l. 

Pitch  and  Lateral  channels  are  coupled  via  constant 
missile  role  rate  P  .   Also,  the  linear  model  has  inertial, 
kinematic  and  aerodynamic  cross  couplings,  which  will  be 
discussed  further,  in  the  non-linear  model. 

1.  Pitch  Channel 

The  linear  aerodynamic  model  of  the  pitch 
channel  consisted  of  the  uncoupled  model  mentioned  in 
chapter  II,  which  is  coupled  with  the  Roll-Yaw  channels  via 
constant  roll  rate  (P  ) ,  as  shown  in  Figure  3.1. 

In  this  work,  the  Pe  was  set  to  zero,  because 
when  Pe=0,  the  Qe  (i.e  equilibrium  pitch  rate)  has  been 
found  to  have  negligible  influence  in  the  lateral  model 
(i.e  roll-yaw) . 

The  equations,  which  represent  the  pitch 
aerodynamic  model  given  in  [Ref.  2]  are  the  same  as  the 
equations  of  the  uncoupled  linear  model,  mentioned  in 
chapter  II. 

2.  Lateral  Channels 

The  linear  lateral  aerodynamic  model  consisted 
of  the  uncoupled  roll  and  yaw  aerodynamic  channel,  coupled 
via  the  linearized  aerodynamic  derivatives  C    ,  C    and 

C^   .  Also,  it  is  coupled  with  the  pitch  dynamic  model  via 
5R 
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the   Pe,    as   previously  mentioned.      The    lateral    dynamic   model 
is    shown    in   Figure    3.2,    as   given    in    [Ref.    2]. 

To    find    the    equations   which    represent    the 
lateral   dynamic   model,    Pe    and    Qe    are    set    to    zero    because    of 
their    negligible    influence    in    the    lateral    channel. 

The    aforementioned    equations    are    as    follows: 
*Wz    q-S-£<  /w     (8««)  (II. B. 2-1) 

E<=  6.CY4-hCyjy.<5y  Ut%)  (III. B.  2-2) 


where : 


t  rSWqi-J-Es/l**  (deglsec2)  (iii.b.2-3) 

h-  tf-Citf  4  GnJp  .£*+  GuJy  Jy  (d«$  (III  .B.  2-4) 

pr  5^3^^d  £q /Ixx    Uegfsec2)  (in.B.2-5) 

4>-  P/5T.3      CrotdlseO  (iii.b.2-6) 

E3T  Cyy  -<Jy  +  Qfg*<?p   Cdeg)  (iii.b.2-7) 

E14 ■=■  £3f  Ce^.S      (deg)  (in. b. 2-8) 

6-  BA|5f.3-r+\<.(q^M  (deglsct)  (in. b. 2-9) 


where: 

k 

= 

0.48 

S 

= 

*(ft)2 

d 

- 

2.0(ft) 

q 

= 

1650    (lbf/ft2) 

w 

= 

2525    (lbf) 

Ixx      =    40.0    (slug-ft2) 
Izz      =    810.0    (slug-ft2) 

cY 

6Y   =      0.22 


91 


IT) 

<u 
Q_ 

a) 

TJ 

o 

s 

u 

•r-1 

s 

(TJ 

c 

>1 

T> 

, , 

? 

<TJ 

5H 

\ 

r-H 

r-( 

. 

O 

TJ> 

, . 

cc 

V 

u 

s./ 

^— 

T3 

<U 

o 

CO 

1) 

J£ 

o> 

r-H 

*"*^ 

(J 

cu 

<U 

o> 

ai 

■a 

S_i 

<u 

TJ 

*-> 

1-" 

(1) 

«J 

■P 

<U 

— 

TO 

as 

O 

0) 

JZ 

r— 1 

a 

zz 

O) 

+-> 

-. 

O 

c 

(0 

a 

(0 

*-) 

<-• 

*-> 

1) 

c 

c 

c 

c 

<T3 

"5 

<o 

•■H 

♦J 

in 

.J 

C 

c 

C 

o 

o 

o 

u 

o 

u 

II 

CN 

II 

II 

• 

a. 

<? 

0) 

a 

• 

' 

CM 

CO 

0^ 

•r-t 

0) 

CL4 

+-> 

O 

92 


cY 
xe  =  -0.082 


cn 

nSR  =   0.018 

c 

n6   =  -0.019 


C 

6y  =  -0.053 


Cn 


C* 


6Y  =  -0.016 


3   =  -0.009 

6R  =  0.035 
A     =  10° 
As  referred  to  in  Appendix  A  and  Table  I. 

Substituting  the  values  into  the  above  equations, 
they  yield: 

T-  -13.93M-S+ 13.95*1  Sp-  3S.S7-&Y (d^trfj (iii.b.2-10) 
p  I  -  133.66 •*+519.T?8-^-33T-Wiy(*8b^)(III-B- 2-lD 
g  rO.onM5P-T-O.0S0S84.O.aif5v(deg/sec)(III-B-2"12) 
4>  =    P/5T.3    (radlseO  (in. b. 2-13) 

C.   CBTT  AUTOPILOT  CONTROL  LAWS  FOR  CIRCULAR  AIRFRAME 

The  control  laws,  which  were  used  by  the  CBTT 

autopilots  of  the  circular  and  elliptical  airframe,  are 

shown  in  Figures  2.6  and  3.3,  as  given  in  [Ref.  2]. 

The  pitch  control  law  for  the  circular  airframe  is  the 

same  as  determined  in  the  uncoupled  pitch  channel  study 
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section  II. B. 2,  and  therefore,  the  equations  which 
represent  it,  are  the  same  as  in  section  II. B. 3.,  except 
the  equation  for  pitch  tail  incidence  angle  which  is: 

J"Pz  -1g&4-<£p-M<H95.3-<$pc  (devisee)    (in.c.i) 

The  uncoupled  roll  and  yaw  laws  are  coupled  via  a 
cross-coupling  branch,  shown  in  solid  line  in  Figure  3.3, 
in  order  to  consist  the  linear  coupled  lateral  control  law, 

In  the  roll  actuator  of  the  lateral  control  law,  the 
frequency  has  been  changed  from  110  rad/sec  to  60  rad/sec, 
for  improving  of  autopilot  stability. 

The  autopilot  cross-coupling  branch  between  Roll  and 
Yaw  channels  has  been  added  to  provide  coordinated  motion. 
Coordinated  motion  or  zero  sideslip  angle  is  achieved  by 
directing  the  body  fixed  pitch  axis  of  the  missile  of  the 
missile  velocity  vector,  so  that  there  is  no  component  of 
missile  velocity  along  the  body  fixed  yaw  axis  of  the 
missile.   Figure  3.4,  as  given  in  [Ref.  2]  shows  the 
attitude  of  the  missile  body  with  respect  to  its  velocity 
vector  V. 

When  commanding  an  upward  maneuver  (i.e,  *  =0),  the 
missile  body  moves  upward  with  its  pitch  axis  directed  at 
the  velocity  vector  until  it  reaches  the  desired  maneuver 
level  or  angle-of-attack.   No  roll  motion  is  required  to 
maintain  coordination  for  this  maneuver. 
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Fig.  3.4   Coordinated  Missile  Motion  for 
Coordinated  bank-to-turn  (CBTT) 
and  limited  bank-to-turn  (LBTT) 
Control  policies 
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For  maneuvers  in  the  *=45  or  90  degree  directions,  the 
missile  body  rolls  the  desired  *  and  moves  its  pitch  axis 
directed  at  the  velocity  vector,  until  it  achieves  the 
desired  maneuver  level.   For  maneuvers  with  *  higher  than 
90°,  the  missile  must  roll  about  the  velocity  vector,  while 
the  Yaw  channel  directs  the  pitch  axis  towards  the  velocity 
vector  for  minimum  sideslip  angle. 

The  missile  is  forced  to  roll  about  its  velocity 
vector,  when  the  desired  maneuver  direction  is  in  the 
negative  angle-of-attack  direction,  because  it  is  desired 
to  avoid  negative  angles-of-attack.   These  maneuvers  are 
illustrated  in  Figure  3.4. 

The  coordinating  command  r  ,  is  a  yaw  angular  rate  from 
a  rate  gyro.   The  command  is  equal  to  -  Pa_,  where  a  =a 
(Figure  3.3)  in  the  coordinating  branch  is  the  equivalent 
angle-of-attack,  which  is  exactly  equal  to  angle  of  attack 
(a) ,  in  the  linear  studies.   The  equations  which  represent 
the  lateral  channel  are  as  follows: 

1.   Roll  Channel 

a.  Roll  Command  Compensator  Equation 
Same  equation  as  in  section  II. E 

Xr  ?X  -HTG  (<J>c-4>)  (in. c. l-i) 

b.  Roll  Rate  Compensator  Equation 

Si  =  — : (x-  p/s*3) 

1      S    ,  '     (III. C. 1-2) 
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Inverse  Laplace  transforming  and  rearranging 
the  equation  (III. C. 1-2),  it  becomes: 

Y,  =  -5t,  +  0. 0503  (X-P/5T-3)  t  0.*55  (X-  P/5f.3)    (III. C.  1-3) 

c.  Psuedo  Differentiator  Equation 

X. 0-015- S (  p/5T3^       (III. C. 1-4) 

-1T-  +  4 

Utilizing  inverse  Laplace  transformation  and 
minor  manipulating  the  above  equation,  one  has: 

Xt  r  -6X1  +  0.072  (  Pi  513)  (III. C.  1-5) 

d.  Actuator's  Compensator  Equation 

The  frequency  of  the  compensator  was  changed 
from  110  rad/sec  to  60  rad/sec  for  purposes  of  stability, 
as  referred  in  section  III.C.   The  equation  represents  the 
network  is: 

5     +i 


$Rc- — (y-x/) 


60 

(III.C. 1-6) 

15 
Using  inverse  Laplace  transformation,  the 

equation  III.C. 1-6,  one  obtains: 

So   =   -15  SRc  +  0-25  (Y-X,)  +15  (Y-X,)  (iii.c. i-7) 

e.   Actuator  Equation 

fRr-IIi_  <fcc  (^        (iii.ci-8) 

Utilizing  inverse  Laplace  transformation  and 
rearranging  the  equation  III.C. 1-8,  it  yields. 
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t 


(III. C. 1-9) 

All  the  referred  to  above  state  variables,  are 
shown  in  figure  3.3. 
2.   Yaw  Channel 

a.  Acceleration  Compensator  Equation 

\f  .     0.32 ^  (III. C. 2-1) 

Inverse  Laplace  transforming,  equation 
III. C 2-1  becomes: 

Y  -  -5Y  +1.6  Tly  (III. C. 2-2) 

b.  Actuator  Compensator  Equation 

5ycr  (<tf5(^-M)/$).  (Y-0.M58a.p/5T.5+t)   (III. C. 2-3) 
Utilizing  inverse  Laplace  transformation   and 
rearranging  the  equation  III. C. 2-3,  it  yields: 

<$YC  -  O.^5(Y+y(5T.3-0.M5Sa.F>/5r3)+  M.85(y+r|5?.3-a«8.a.F/5T.3^ 

( I II. C. 2-4) 

c.  Actuator  Equation 

<ty= -188.4 -Jy  -H0?9S.3cVc  (III. C 2-5) 

Making  use  the  state-variable  representation, 
the  above  equations,  of  the  aerodynamic  models  and  control 
laws,  can  be  modelled  in  the  following,  state  variable 
equation  system: 
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(III. C. 2-6) 


where    the    state   variables    are: 


r 
p 

3 
X 

Y- 


yaw  angular  rate 

roll  angular  rate 

roll  angle 

sideslip  angle 

output  of  roll  command  compensator 

output  of  roll  rate  compensator 

output  of  pseudo  differentiator 

input  command  in  roll  actuator 

roll  tail  incidence 

output  of  yaw  acceleration  compensator 

input  command  in  yaw  actuator 

yaw  tail  incidence  angle 

All  the  referred  before  state  variables  are 
shown  in  Figures  3.2  and  3.3. 
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R 


Y 

6, 
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D.   ANALYSIS  OF  LINEAR  CBTT  AUTOPILOTS  FOR  CIRCULAR 

AIRFRAME 

In  the  analysis  of  the  linear  CBTT  autopilots  it  is 
assumed  that  the  missile  is  initially  in  the  desired 
maneuver  plane,  is  trimmed  at  ten  degrees  angles-of-attack, 
the  equilibrium  roll  rate  Pe  and  equilibrium  pitch  rate  Qe 
were  set  equal  to  zero. 

For  purposes  of  analysis,  a  CSMP  program  was  written 
(Appendix  G)  using  the  equations  of  the  aerodynamic  models 
and  control  laws  of  the  linear  CBTT  autopilots  for  the 
circular  airframe. 

Figure  3.5  shows  the  sideslip  angle  for  the  linearized 
CBTT  autopilot,  having  a  maximum  value  of  2.32°.   Since  the 
model  is  linear,  and  therefore  the  magnitude  of  the 
sideslip  is  directly  proportional  to  the  magnitude  of  the 
input  command,  a  one  radian  roll  angle  command  was  used  for 
convenience . 

The  coupled  autopilots  have  a  roll  angle  response  which 
is  the  same,  as  the  uncoupled  roll  channel  response, 
differing  slightly  in  the  overshoot,  as  shown  in  Figure 
3.6. 

The  aforementioned  values  and  plots  are  exactly  the 
same  as  the  plots  given  in  the  [Ref.  2]  therefore  its 
results  are  verified. 
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Fig.  3.5   Sideslip  angle  (3)  vs  Time  (t) 

Linear  CBTT  Autopilot  (1  RAD  Roll 
Command);  Circular  airframe  (a  =10  deg , 
Pe=Qe=0  deg) 
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Linear  CBTT  Autopilot  (1  RAD  Roll 
Command)  Circular  airframe  (a  =10 
deg ,  Pe=Q  =0  deg) 
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Comparing  the  state  variable  form  system  with  these 
previous  discussed  in  the  section  II. D  and  II. E,  this  is 
the  twelfth  order  system  because  it  includes  roll  and  yaw 
autopilots  and  the  matrices  A,B  are  different  than  the  A,B 
matrices  at  the  uncoupled  roll  and  yaw  autopilot,  because 
of  the  existance  of  the  coordination  branch. 

E.   AERODYNAMIC  MODELS  FOR  ELLIPTICAL  AIRFRAME 

The  dynamic  models  for  elliptical  airframe  are  the  same 
as  those  of  circular  airframe  and  are  shown  in  Figures  3.1 
and  3.2. 

The  differences  of  the  circular  and  elliptical  models 
are  in  the  following  constants  and  linearized  aerodynamics 
derivatives: 

w     =  2475  (lbf) 

Ixx   =  110  (slug-ft2) 

Izz   =  790  (slug-ft2) 

CY 

«Y   =   0.015 

CY 

e   =  -0.054 
c 

n5D   =   0.014 
K 

c 

8    =   0.024 

C 

6y   =  -0.039 
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6y       =    -0.010 

a       =  -0.027 

6R      =    0.023 

as    referred    to    in   Appendix   A   and   Table    I.      The    equation 
which    represent    the   aerodynamic   model,    are   as    follows: 

f    Z    1S.0M69  -  6  +-I0.52qr3-<5"R  -  29.32^2-<?V  (in.E.l) 

p    r  --H5.W)*.S+  12q.}09<-<$B-5M.oq.Sy  (iii.E.2) 

4    -     P|573  (III.E.3) 

$Z    0.HM52  p-r  +  0.05UWt?-O.Of5W<fy         (III.E.4) 

F.   CBTT  AUTOPILOTS  CONTROL  LAWS  FOR  ELLIPTICAL  AIRFRAME 

The  control  laws,  which  were  used  by  the  CBTT 
autopilots  of  the  elliptical  airframe  are  the  same  as  those 
of  the  circular  airframe  and  are  shown  in  Figure  2.6  and 
3.3. 

The  differences  which  exist  in  the  control  laws  of  the 
two  airframes  are  the  gain  Kyp  of  the  coordinated  branch, 
and  the  gains  of  the  actuators. 

Hence,  the  equations  that  represent  the  control  laws 
are : 

1.   Pitch  Channel 

The  linear  pitch  control  law  equations  for 

the  elliptical  airframe,  are  the  equations  of  the  Pitch 

uncoupled  channel,  as  given  in  section  II. C. 5. 
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2.   Roll  Channel 

Same  equations  as  in  par.  III.C. 

a.-X  Z-8X+H.6  (4>c-4>)  (II. C. 1-1) 

b.  ^--5^+0.0505  (x-P|5T.3)+0.?55U-P/57.3)  (II. C. 1-3) 

c  X^r  -6X,  +  0.0?8  (P/5r.3)  (II. C. 3-5) 

d.      Actuator    Compensator 

(III. F. 2-1) 


S^JfcwJAlil  ( w\ 


Utilizing    inverse    Laplace    transformation   and 
rearranging    the   equation   III. F. 2-1,    it   becomes: 

Src  =  -15^Rc  + 10^25 (VX|)  "+  62.53  (V*<)  (III. F.  2-2) 

e.     5rz  -ig?.q£fc  +  10.?<J53;c1o3Src  (III.C. 3-8) 

3.      Yaw  Channel 

a.  Acceleration  Compensator  Equation 

w            0.839         ^  (III. F. 3-1) 

y  =— vi  y 

Inverse    Laplace    transforming   and    rearranging 
the   equation   III. F. 3-1,    it   yields: 

Y:  -*IV  +3.556  Hy  (III.F.3-2) 

b.  Actuator  Compensator  Equation 

JYc  _  0.6Ot,l+  i.OS        f  Y  +  Wsf  3  _  -  p/5T.5) 

(III. F. 3-3) 
Utiliing  inverse  Laplace  transformation  and 
rearranging  the  above  equation,  one  obtains: 
<$YC  =  0.608  (Y4^-2£)  +6-08  (y^-JL  )    (III.P.3-4) 
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c.   Actuator  Equation 

Same  as  in  par.  III.F.2.e 

Utilizing  the  state-variable  representation, 
the  equation  which  represent  the  aerodynamic  model  and  tne 
control  laws,  can  be  modelled  in  the  following 
twelfth-order  system: 
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where  the  state  variables  are: 


8 
X 
Y. 


yaw  angular  rate 

roll  angular  rate 

roll  angle 

side  slip  angle 

output  of  roll  command  compensator 

output  of  roll  rate  compensator 

output  or  psuedo-dif f erentiator  network 

input  command  in  roll  actuator 


<*c 


R 


roll  tail  incidence  angle 
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Y   :     output  of  acceleration  compensator  in  yaw 
control  law 

S   :     input  command  in  yaw  actuator 
c 

5   :     yaw  tail  incidence  angle 

All  the  state  variables  ace  shown  in  Figure  3.2  and 
3.3.   The  above  system  is  different  in  the  matrix  A  and  B 
from  the  system  discussed  in  chapter  II  for  uncoupled  Roll 
and  Yaw  autopilots,  because  of  the  coordination  branch 
equation . 

G.   ANALYSIS  OF  LINEAR  CBTT  AUTOPILOTS  FOR  ELLIPTICAL 

AIRFRAME 

In  the  analysis  of  the  linear  CBTT  for  elliptical 
airframe,  the  same  assumptions,  which  were  made  for  the 
circular  airframe,  are  valid. 

A  CSMP  program  was  written  (Appendix  G)  with  the  above 
equations,  which  represent  the  aerodynamic  models  and  the 
control  laws  of  the  linear  CBTT  autopilots. 

Figure  3.7  shows  the  sideslip  angle  8  for  the 
linearized  CBTT  autopilot,  having  maximum  values  of  -0.46 
and  0.26°. 

Figure  3.8  shows  the  roll  angle  response  which  has  an 

overshoot  of  4  percent  and  is  about  the  same  with  the  roll 

angle  of  uncoupled  roll  autopilot.   The  above  plots  verify 

the  results,  given  in  the  [Ref.  2]. 
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Fig.  3.7   Sideslip  angle  (B)  Ys  Time  (t) 

Linear  CBTT  Autopilot  (1  RAD  Roll 
command)  Elliptical  airframe  (a   =10°, 
Pe=Qe=0°) 
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H.   RELATIVE  STABILITY  OF  CBTT  AUTOPILOTS  OF  CIRCULAR  AND 

ELLIPTICAL  AIRFRAMES 

When  the  coordinating  branch  gain  Kyp  was  set  to  zero 
in  CBTT  autopilot  of  the  circular  airframe,  the  maximum 
sideslip  angle  increased  to  2.8  degrees.   Therefore,  for 
the  circular  airframe  the  coordinating  branch  is  not  very 
effective  in  helping  the  yaw  autopilot  to  reduce  sideslip 
angle . 

If  kyp  is  increased  to  unity,  as  it  is  for  the  CBTT 

autopilot  of  the  elliptical  airframe,  the  maximum  sideslip 

remains  below  one  degree  for  3  seconds  but  the  autopilot  is 

unstable.   Modifying  of  the  circular  airframe  will  be  done 

in  the  non-linear  CBTT,  for  minimization  of  the  sideslip 

angle . 

Removal  of  the  aerodynamic  cross-coupling  (i.e  C    , 

B 
C    ,   C    )  had  no  effect  on  the  sideslip  angle.   This 

Y       6R 
showed  that  the  aerodynamic  cross-coupling  plays  an 

indirect  role  in  determining  the  quality  of  sideslip 

control  by  determining  the  relative  stability  of  the 

coordination  branch  or  the  magnitude  of  the  coordination 

gain  K   .   The  magnitude  of  sideslip  angle  is  dependent  on 

the  nulling  effects  of  two  parallel  paths  shown  in  Figure 

3.9,  as  given  in  [Ref.  2]. 
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The  contribution  of  yaw  acceleration  n   to  the  maximum 
0  is  negligible  and  therefore  neglected.   3  is  formed 
mainly  by  the  subtraction  of  the  kinematic  paths  of  a  p  and 
the  yaw  angular  rate  r. 

The  coordination  is  obtained  in  the  CBTT  control  law  by 
commanding  the  yaw  autopilot  with  a  yaw  angular  rate 
command  rc  (Figure  3.2)  of  Kyp  aep,  which  forces  r  to  be 
equal  to  ae.p  and  therefore  nulling   as  shown  in  Figure 
3.9.   The  nulling  process  will  be  accomplished  more 
efficiently  if  Kyp  =1.0.   The  reason  the  sideslip  is  not 
nulled  completely  is  that  n   is  not  zero  and  r  can  not 
equal  r   over  all  frequencies. 

The  linearized  CBTT  autopilot  of  the  circular  airframe 
was  unstable,  when  the  uncoupled  channel  control  laws  used 
with  a  coordinating  branch  gain   KYp  of  unity.   The 
autopilot  was  stabilized  by  decreasing   to  0.458.   The 
actuator  of  roll  control  channel  has  a  frequency  of  60 
rad/sec. 

The  stability  of  the  CBTT  autopilot  for  the  elliptical 
airframe  is  considerably  better,  due  to  the  setting  of  KYp 
at  unity,  which  results  in  a  substantial  decrease  in  the 
magnitude  of  sideslip  angle  as  shown  in  Figure  3.9. 
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I.   CONCLUSIONS 

1.  The  linear  coupled  autopilots  verify  the  results  of 
the  uncoupled  autopilots  and  keep  the  sideslip  angle  small 
for  coordination  motion,  (section  II. D.  and  II. E) 

2.  Further  Conclusions  about  the  relative  stability  of 
the  autopilots,  about  the  gain   Kvp  of  the  coordination 
branch  and  about  the  sensitivity  of  the  autopilots  to  the 
aerodynamic  cross-coupling  are  given  in  chapter  7  of  [Ref. 
2]. 


114 


IV.   NONLINEAR  ANALYSIS  OF  CBTT  AUTOPILOT  FOR  CIRCULAR  AND 
ELLIPTICAL  AIRFRAME 

A.    GENERAL 

The  three  dimensional  nonlinear  aerodynamic  models 
used  for  the  analysis  are  shown  in  the  Figures  4.1  and  4.2. 
For  the  configuration  under  consideration,  the  assumptions, 
which  have  been  found  to  be  consistent  are  referred  to  in 
chapter  5  of  [Ref.  2], 

Nonlinear  aerodynamics  and  mass  parameters  values  are 
presented  in  Appendixes  A  and  B,  as  taken  from  the  [Ref. 
2].   The  same  flight  conditions  (i.e  60kft  altitude  and 
M=3.95)  as  in  linear  CBTT. 

In  this  section  the  control  laws  for  the  non-linear 
CBTT  autopilots  are  discussed.   The  results  for  commanding 
the  CBTT  autopilots  for  both  airframes  are  shown  and  the 
desired  aerodynamic  model  to  enhance  CBTT  performance  is 
determined.   The  performance  of  the  CBTT  autopilot  for  the 
elliptical  airframe  is  determined,  when  no  lateral 
aerodynamic  cross-coupling  exists  and  with  ideal  airframe 
dynamics.   Also  the  performance  of  the  ?qCBTT  autopilot  for 
the  circular  airframe  is  determined,  when  no  inertial  and 
kinematic  cross-coupling  exist  into  the  pitch  autopilot. 
The  conclusion  of  the  nonlinear  analysis  are  stated. 
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B.    CONTROL  LAWS  FOR  ELLIPTICAL  AIRFRAME 

The  control  laws  used  for  the  following  nonlinear  3-D 
studies  were  the  same  as  those  used  for  the  linear  studies 
in  section  III,  which  are  shown  in  Figures  4.3  and  4.4 
taken  from  [Ref .  2] ,  except  for  the  gain  as  shown  in  the 
bold  line  of  the  coordination  branch  in  Figure  4.4.   The 
new  gain  a  is  held  constant  at  one  degree  magnitude  for 
angles  of  attack  less  than  one  degree  positive  and  greater 
than  negative  5  degrees.   For  angles-of-attack  greater  than 
one  degree  positive,  the  gain  a  is  equal  to  the  angle  of 
attack.   This  maintains  coordination  for  very  small 
angles-of-attack. 

Gravity  effects  were  not  included  in  the  linear 
studies,  because  it  was  assumed  to  have  a  negligible 
influence  on  autopilot  stability  and  response  perturbations 
about  a  missile  trim  conditon. 

Gravity  effects  were  included  in  the  nonlinear 
studies,  where  the  missile  body-fixed  yaw  axis  will  be 
subjected  to  the  full  force  of  gravity  and  may  therefore 
have  a  significant  influence  on  sideslip. 

1.   Pitch  Control  Law 

The  pitch  control  law  for  the  elliptical  airframe 
is  the  same  that  is  used  in  the  linear  studies  in  section 
III.F.,  except  in  the  command  applied  to  the  channel. 
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In  inertial  rectangular  coordinates,  the 
acceleration  commands,  are: 

a.  n„  (acceleration  command  in  inertial  Zv 

c  v  Y 

direction  as  in  figure  2.4)  =  nz   -  cos*    IV. B. 1-1 

c 
where:   n7   :    guidance  command  (gees) 
c 

-cos*  :    anti-gravity  bias  command  (gees) 

9   :    Elevation  Euler  angle  =   (qcos^-rsin*) dt 

b.  ny   acceleration  command  in  inertial  Y 

c 

direction,  which  is  equal  to  zero,  because  it  is  desired 
coordinated  motion.   In  the  polar  coordinates,  the 
acceleration  commands  are: 

a.  <vic  r  -  [  (^2C-  coae)*^2  ]  z 

z  -  (i\zc-  cose) 


b-     <K   _  L-J    /*N 


*c=W(^) 


(IV. B. 1-2) 


(IV. B. 1-3) 


hence  *   =  nn  and  n  =0,1,2... 

Since  the  pitch  control  law  of  the  elliptical 

airframe  does  not  have  an  integration  in  the  acceleration 

error  path  (i.e  in  the  acceleration  compensator  as  in  the 

case  of  circular   airframe)  and  requires  a  gain  in  series 

with  the  acceleration  command  shown  in  Figures  4.3  it  was 

necessary  to  modify  the  antigravity  command  as  follows  to 

assure  an  anti-gravity  bias  of  just  one  gee. 
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n     (elliptical  airframe  inertial  acceleration  command)  - 

c        c 

-  n  -  0.913  cose  (III. B. 1.4) 

c 
The  autopilot  was  tested,  using  a  commanded  2  gee 

(0°,  180°),  that  denotes  an  inertial  command  of  2  gees, 

which  is  first  applied  in  the  0   or  upward  direction  at  t  = 

2  seconds  for  3  seconds  and  at  the  fifth  second  is  applied 

in  the  180   or  downward  direction.   Since  at  first  both  the 

missile  roll  angle  and  roll  angle  command  are  at  zero 

degrees,  there  is  no  roll  motion  and  the  missile  turns 

upward  as  a  skid-to-turn  controlled  missile.   At  fifth 

second  the  missile  is  commanded  to  roll  through  180  degrees 

while  moving  in  a  coordinated  manner  in  yaw  and  roll  to 

minimize  sideslip  angle  and  prevent  or  minimize  negative 

angle-of-attack. 

The  other  equations,  which  represent  the  nonlinear 
pitch  control  law  are  the  same  as  for  linear  studies 
section  III.C.  and  the  state  variables  are  shown  in  Figure 
4.3. 

The  equations  are: 

X   r -150X -M50%B  (iv. b. l-i) 

<f    r  -6Y-  0.5305  1\c  +  OMX  (iv.B.l-2) 

<Jpc  -  -0.383*  (Y+Q|5?.5)-3.0*(Y-9l5?.3)        (lV.B.1-3) 

<5f  =  -AS8.H    Sp  + 40^5.3 -&pc  (iv. b.  1-4) 
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where : 


n7    :   achieved  body-fixed  acceleration  in  Z  axis 

n„   :    inertial  acceleration  command, 
c 


2.   Lateral  Control  Laws 

The  lateral  control  laws  are  the  same  as  those 

used  for  the  linear  studies,  except  the  modification  in  the 

coordination  branch. 

The  command  $   for  the  roll  control  law  is  zero 
c 

for  the  first  5  seconds  and  180   affterwards,  as  derived  in 
IV.B.l. 

The  equations  which  represent  the  lateral  control 
laws  are  the  same  those  of  section  III.F. 


Y2  r  -q  Y2 -1-3.556  >ty6  (iv. b. 2-1) 

Si,z  0.6O8  (v4-i^')  +  ^*ry, +  — -^1 

c  2       5T.3  .   (57-3)*  '  V   2        57-3       C57-3)2  / 

V2    57.3    5X3.23     3233.23  r6'0SV2  +  ^-  3233^) 


(IV. B. 2-2) 


where: 


and 


a>l° 


dc-  d  J 

<M*°l  a<i° 
acr  0  J 

l,  =  -8JCl»1U*4.i«*c  (IV.B.2-3) 

<£y  ~  -lMM  Sy  +  i0}95.3  <$YC  (iv.b.2-4) 

%  z  -51ff +0.05033(irP/57.3)  +  0J55  (X«-Pfe*3) 

X2-  -6Xa  + 0-0*3  (P/5«} 

S^r  -15  &c  4iP425CY4-XO+62-55  6fl-Xi)      (iv.b.2-7) 

c  123 


(IV. B. 2-5) 
(IV. B. 2-6) 


The  state  variables  are  shown  in  Figure  4.4.  For 
purposes  of  analysis,  a  CSMP  program  was  written  (Appendix 
H). 

C.    AERODYNAMIC  MODELS  FOR  ELLIPTICAL  AIRFRAME 

The  aerodynamic  models  of  Pitch  and  Roll-Yaw  channels 
for  the  elliptical  airframe  are  shown  in  Figures  4.1,4.2. 
The  airframe  dynamics  under  consideration,  have  a  level  of 
complexity  sufficient  to  determine  the  critical  areas  of 
concern  regarding  the  stability  and  response  of  CBTT 
control . 

1.   Pitch  Channel 

The  nonlinear  aerodynamic  model  of  the  Pitch 
channel  consisted  of  the  uncoupled  model,  coupled  with  the 
Roll-Yaw  channels  via  the  inertial  cross-coupling  of  rp 
into  q  and  the  kinematic  cross-coupling  of  -B.P  into  a. 
Also  the  antigravity  command  cos*  cose  is  inserted  into  n  . 
This  model  is  shown  in  Figure  4.1. 

The  above  terms  consist  of  the  nonlinear  terms  of 
the  pitch  channel. 

The  equations  which  represent  the  pitch 
aerodynamic  model  are: 

%  r  -q^CN/w  (g*es)  (iv.c.i-i) 

q   r  p.r/5*.3+5*.3q£dCMllyy  (iv.c.i-2) 
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d  nq+K.^-  ?.e/5ll  (iv. c. 1-3) 

where : 

1845 

K    =  =  0.48 

V 

q   =  1650  (lbf/ft2) 

S    =  ir  (ft2) 

d    =2  (ft) 

w    =  2475  (lbf) 

Iyy  =  790  (slug-ft2 

Functions  CN(a,s  )  and  cm(a'6D)  are  nonlinear 
functions  shown  in  Appendix  B  that  vary  with  a  and   6p. 

Substituting  the  values  into  the  above  equations, 
they  yield: 

^=  -J.3333-C,n/  (IV.C.1-4) 

q  =  0.mS2-M  +  1SL9S4J   Cm  (iv. c. 1-5) 

i  -  q-0.639<3*Ctf-0-CMq52p.€       (iv.c.i-6) 

2.   Lateral  Channels 

The  nonlinear  aerodynamic  model  consisted  of  the 

uncoupled  roll  and  yaw  channels,  which  coupled  via  C 

B 
C    and  C_   .   Also  the  lateral  channel  coupled  with 

the  lateral  channel  coupled  with  the  pitch  channel  via  the 
inertial  cross-coupling  of  -qp  into  r  and  the  kinematic 
cross-coupling  of  ap  into  B.   The  above  term  is  the 
nonlinear  part  the  lateral  pitch  channel.   The  nonlinear 
lateral  dynamic  model  is  shown  in  Figure  4.2.   The 


125 


(IV.C.2-4a) 
(IV. C. 2-5) 


equations    that    represent   the   aforementioned   model    are   as 
follows : 

V  =  S^  (<*«•*  +  C^y  •**)  /  W     <>es)  (iv. c. 2-1) 

6    =  K*WB  +  a.pJsfS  -r    (<te#lsec)  (iv.c.2-2) 

(IV. C. 2-3) 

P    =   5T.5  9^<i    (<««•£  +   Cejy  •J/+Qfc.&)  /ixx       (iv. c. 2-4) 
<i>    -     P/5T-3    (rouJUec) 
4  =    P     (dtglsec) 

where: 

K  =  0.48 

q  =  1650    (lbf/ft2) 

S  =  it    (ft2) 

d  =2    (ft) 

w  =  2475    (lbf) 

Ixx  =  110    (slug-ft2 

I        =   853    (slug-ft2) 
zz 

as   referred    to    in  Appendix  A. 

Functions   Cy       (a) ,    Cy    (a) ,      Cn    (a) ,      Cn       (a) ,    C£ 

5y       B        B        6y       «y 

Cj.   ^  '  c«.   ^a^  and  c    ^  are  nonlinear  functions  which 

8        SR  SR 

vary  with  a  angle-of-attack  and  referred  in  Appendix  B. 

Substituting  the  values  into  the  referred  before 

equations  one  obtains: 
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^bz  1.3333  (c1€-8  +Cyjy.5y) 
6  =  0.63998  (%-8  ±0itrfo)+Q.0mSlCL'P-Y 

p  -  5M00.3OT  (C^.S+C^-fy  +  C^.^) 


(IV. C. 2-5) 
(IV. C. 2-6) 
(IV. C. 2-7) 
(IV. C. 2-8) 
(IV. C. 2-9) 


D.    ANALYSIS  OF  NONLINEAR  CBTT  AUTOPILOTS  FOR  ELLIPTICAL 

AIRFRAME 

For  purposes  of  analysis,  a  CSrtP  program  was  written 

(Appendix  H)  with  the  aforementioned  equations  of  the 

aerodynamics  models  and  control  laws  of  the  nonlinear  CBTT 

autopilots. 

To  decrease  computer  time,  appropriate  initial 

conditions  were  inserted.   The  initial  conditions  are: 

a   :  angles  of  attack  =  2.41  degrees 

5   :  pitch  tail  incidence  =  0.658  degrees 

9   :  pitch  Euler  angle  =  0.636  rad 

X   :  output  of  pitch  acceleration  feedback  lag  =  lgee 

Y   :  output  of  pitch  acceleration  compensator  =  -0.0105 

5   :  pitch  actuator  command  =  0.0148  rad 
P 

The  achieved  maneuver  plane  acceleration  n   is 

z 

calculated    from   the   body-fixed    accelerations   n_    (Figure 

8 
4.6)    and      (Figure    4.7)    and    the    roll    angle    <j>    (Figure    4.13) 

as  follows: 
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Fig.  4.5  Achieved  Inertial  Acceleration  (nz)  vs  Time  ^t) 
CBTT  of  Elliptical  airframe:  2  Gees  (0  ,  180  ) 
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Fig.  4.6   Achieved  Body-fixed  Acceleration  (n„  )  vs 

B 
Time  ( t)  CBTT  of  Elliptical  airframe;  2  Gees 
(0°,180°) 
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Fig.  4.7  Achieved  Body-fixed  Acceleration  (ny  )  vs 


B 


Time  ( t)  CBTT  of  Elliptical  airframe;  2 
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Yl   -  ^    .   C0S4>  4-,nYB*  S%IY\4?  (IV.D-1) 

Figure  4.5  through  4.16  show  the  responses  of   n  , 

nY  /  a/  Q  f    <3#  P'  r'  *»  5D'  6V  and  6R*   Also  these  figures 
and  table  II  show  the  time  constants  and  the  overshoots  of 
the  referred  before  variables. 

Figure  4.5  shows  that  the  achieved  inertial 
acceleration  n   doesn't  satisfy  the  requirement  of  the 
overshoot  (i.e  below  10  percent)  and  it  is  not  acceptable. 
The  results  of  this  analysis  verify  the  result  of  the  [Ref. 
2]. 


E.    CBTT  PERFORMANCE  WITH  IDEAL  AIRFRAME  DYNAMICS   FOR 

ELLIPTICAL  AIRFRAME 

The  purpose  of  the  following  simplifications  to  the 
airframe  dynamics  model  is  to  isolate  the  critical 
cross-coupling  paths,  which  have  caused,  the  transients  in 
the  maneuver  plane  acceleration  responses. 

The  dynamic  model  without  the  coupling  paths  will  be 
referred  to  as  ideal  dynamics.   Although  ideal  dynamics  are 
not  physically  attainable,  it  is  a  useful  goal  for  both 
autopilot  and  airframe  designers. 

The  same  guidance  commands  are  applied  to  the  CBTT 
autopilot  of  the  elliptical  airframe  as  in  par.  IV.B.l 
(i.e,  2  gees  (0°,  180°)  but  with  lateral  aerodynamic 
cross-coupling  removed.   That  means,  the  aerodynamic 
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Fig.  4.14   Pitch  tail  incidence  («  )  vs  Time  (t) 
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Fig.  4.15  Yaw  tail  incidence  (6y)  vs  Time  (t) 
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Fig.  4.16   Roll  tail  incidence  (5R)  vs  Time  (t) 

CBTT  of  Elliptical  airframe;  2  Gees  (0  $    180  ) 
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cross-coupling    was  removed.   The  major  influence  of  the 
aerodynamic  cross-coupling  in  the  elliptical  airframe 
response  has  been  to  decrease  the  overshoot  in  the  achieved 
maneuver  plane  acceleration  resulting  from  the  second 
guidance  command. 

In  addition  to  aerodynamic  cross-coupling  the 
kinematic  cross-coupling  of  -ep  into  a  and  inertial 
cross-coupling  of  pr  into  q  were  removed.   Therefore,  the 
only  cross-coupling  which  exist  in  the  airframe  dynamic 
model  are  the  kinematic  coupling  of  ap  into  B  and  the 
inertial  coupling  of  -qp  into  r.   There  is  also  the 
autopilot  cross-coupling  of  the  coordinating  command  from 
the  Roll  to  Yaw  channel. 

The  equations,  which  represent  the  modified 
aerodynamic  model  by  removing  the  referred  cross-coupling 
are  same  with  those  as  in  paragraphs  IV.B.l,  IV. B. 2,  IV.C.l 
and  IV. C. 2,  with  the  only  difference  in  the  following 
equations: 

1.   Pitch  Aerodynamic  Model 

C}    Z   5T.3^^d  CM/Jyy      CJeg/sec2) 

-    T51.*54i  CM  (iv.E.i-i) 

and  • 

Z    q-0.G59<38  CN   (d^/sec)  (iv.e.i-2) 
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Lateral  Aerodynamic  Model 

(IV. E. 1-3) 


and 


Ixx 


.^^ — (%.h 


(IV. E. 2-2) 

For  purposes  of  analysis,  a  CSMP  program  (Appendix 

H)  was  written,  using  these  equations.   Figures  4.17 

through  4.24  show  the  time  responses  of  the  n7,  n7  ,  a,  0, 

L         ZB 
q,  t,  5  ,  6V,  6D.   The  above  figures  and  table  II  show  the 

p     I     K 

time  constants  and  the  overshoots  of  the  variables. 

Figure  4.17  and  Table  II  show  that  the  achieved 
inertial  maneuver  acceleration  satisfy  the  requirements  for 
the  time  constants  and  overshoot,  when  the  cross-couplings 
which  were  discussed  before  have  been  removed. 

Figure  4.25  (taken  from  [Ref.  2]  shows  the 
critical  feedback  paths  which  couple  the  pitch  and  yaw 
channels  via  missile  roll  rate. 

The  aforementioned  figures  and  plots  verify  the 
results  of  the  [Ref.  2]. 

F.    CONTROL  LAWS  FOR  CIRCULAR  AIRFRAME 

The  control  laws  used  for  the  circular  airframe  are 

the  same  as  those  used  in  the  linear  studies,  except  for 
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Angle  of  Attack  (a)  vs  Time  (t) 
CBTT  of  Elliptical  airframe;  Ideal 
airframe  dynamics;  2  Gees  (0°,  180°) 
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CBTT  of  Elliptical  airframe;  Ideal 
airframe  dynamics;  2  Gees  (0°,  180°) 


145 


A 

IA 

\ 

K" 

O 

/»  " 

:  \ 

«n 

.  \ 

m" 

y 

V 

1 

-N-     ^^ 

^ 

cs  o 

UJ 

o 

/ 

o 

/ 

i 

o 

■ 

<*> 

it 

1 

u 

o 

o 

t  ' 



0.0 


1.0 


2.0 


3.0 


4.0 

IISECJ 


s.o 


8.0 


7.0 


9.0 


Fig.  4.20 


Body  pitch  angular 
CBTT  of  Elliptical 
airframe  dynamics; 
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airframe;  Ideal 
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CBTT  of  Elliptical  airframe  Ideal 
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Fig.  4.24   Roll  tail  incidence  («Y)  vs  Time  (t) 
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Fig.  4.25   Critical  Feedback  loop  for  CBTT 
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TABLE  II.  CBTT  PERFORMANCE  OF  ELLIPTICAL  AIRFRAME 


VARIABLE 

:  VALUE           TIME  (sec) 

CROSS- 
COUPLING 

T,(sec) 
Overshoot% 

T^sec) 
Overshoot% 

nZmin 

-2.934(gees) 

5.061 

All 

0.43  and 

0.41  and 

nZmax 

1.3S1 (gees) 

5.661 

All 

0% 

11.25% 

nZmin 

-2. 934(gees) 

5.062 

ideal 

0.48  and 

0.37  and 

Zmax 

0.937 (gees) 

5.853 

ideal 

0% 

0% 

amin 

1.35(deg) 

6.10 

All 

a 
max 

6.75(deg) 

3.10 

All 

amin 

2.22(deg) 

6.295 

ideal 

a 
max 

6.75(deg) 

3.10 

ideal 

min 

-0.9(deg) 

5.49 

All 

B 
max 

1.49(deg) 

6.05 

All 

B  i 
min 

1.53(deg) 

5.45 

ideal 

B 
max 

1.07 (deg) 

5.91 

ideal 

qmin 

-10.6 (deg/sec) 

5.639 

All 

q 
^max 

9.5 (deg/sec 1 

2.337 

All 

qmin 

-7.8 (deg/sec) 

5.328 

ideal 

q 
^max 

9. 5(deg/sec) 

2.337 

ideal 

P  . 

min 

-14.115(deg/sec) 

6.600 

All 

max 

282.094 (deg/sec) 

5.387 

All 

P  . 

min 

-14.29(deg/sec) 

6.600 

ideal 

P 
max 

281 .928  (deg/sec) 

5.387 

ideal 

r  . 
min 

-2.32(deg/sec) 

5.90 

All 

r 
max 

33. 1 (deg/sec) 

5.32 

All 

rmin 

-0.178(deg/sec) 

7.036 

ideal 

max 

37. 5  (deg/sec) 

5.29 

ideal 

Note: 

All  the  above  referred  values  are  for  command  2(0  ,  180  ) 

-J.:  62  percent  time  constant  of  achieved  maneuver  plane 
acceleration  due  to  the  first  guidance  command. 

Y2!  63  percent  time  constant  of  achieved  maneuver  plane 
acceleration  due  to  the  second  guidance  command. 
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the  gain  (a)  of  the  coordination  branch  of  lateral  control 
laws,  as  mentioned  in  the  study  for  the  elliptical 
airframe. 

The  same  2  gees  (0°,  180°)  guidance  commands  and 
flight  condition  were  used  for  direct  comparison  with  the 
performance  results  of  the  elliptical  airframe. 

However,  to  decrease  simulation  run  time,  the  first 
guidance  command  and  anti-gravity  bias  were  applied  at  zero 
time  with  no  missile  or  autopilot  initial  conditions. 

1.   Pitch  Control  Law 

The  pitch  control  law  for  the  circular  airframe  is 
the  same  as  that  use  in  the  uncoupled  autopilot  shown  in 
Figure  4.26.   The  command  used  in  the  autopilot  is 
different  than  the  one  used  in  uncoupled  channel  and  it  is 
as  follows: 

Utr-rfcc-cose  (iv.B.i-D 

as  referred  to  in  section  IV.B.l 

The  equations  which  represent  the  nonlinear  pitch 

control  law  are  the  same  as  those  in  section  III.C.  of 

the  linear  pitch  control  law,  using  K,  =  -0.0387 

A2 
The  equations  are: 

X  r- -150X4-  150  ^B  (IV.F.1-1) 

%-   ^-13531  (vic-X)  (iv.F.i-2) 

%z-5\  4-6.572-(t)c-X)  (iv. f.  1-3) 

6pc- - 0/143. Spc-0. 31058(y3+q/5?.3)-2.23Oj.(Y3-qte?-3) 

(IV. F. 1-4) 
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Sf    =  -18g.4  c5p  +  10*95.3  'Oft.    (IV. F. 1-5) 
The  state  variables  of  the  above  equation  are 
shown  in  Figure  4.26. 

2.   Lateral  Control  Laws 

The  lateral  control  laws  for  the  circular  airframe 
are  the  same  as  those  used  for  the  linear  studies,  except 
the  modification  in  a   (or  a)  which  was  referred  to  before. 

The  same  command  is  used,  as  in  the  elliptical 
airframe.   The  equations  which  represent  the  lateral 
control  laws  are  the  same  those  of  section  III.C.l  and 
III.C.2.  and  they  are: 

Y2  -  -5  Yg  +  16  ty8  (IV. F. 2-1) 

r  0.425  [Y2+^-  -o.oooi33(a€-P+deP)J  +  ^5Cy24-^-o.oooi39^.p) 

(IV. F. 2-2) 
where : 


Ot-  oc  £ 

r  a>iG 

de  ~  d  i 

and 


<Sy    -  -m-H  £y-»- 10*95.3  <$VC  (iv.f.2-3) 

^  r  -8Xi-17.e4>4-lT.e4>c  (iv. f. 2-4) 

Y<  r  -5  ^  +  0.05033  (X,-P/5T3)-K).?5 5  (Xf-P/57<3) 

(IV. F. 2-5) 
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X2  -  "6X2  +  0.0¥S  (  P/5?.^  (IV. F. 2-6) 

<&c-  H5&c  + 0.25  (VrX2)-^5(VrX2>)  (IV.F.2-7) 

5fc  =  -iS^J^  -H  0395.3  •£&.  (iv.  f.  2-8) 

The   above    state   variables    are    shown  in   Figure 


4.27. 


G.    AERODYNAMIC  MODELS  FOR  CIRCULAR  AIRFRAME 

The  aerodynamic  models  of  Pitch  and  Roll-Yaw  channels 
foe  the  circular  airframe  are  the  same  as  those  of  the 
elliptical  airframe  and  are  shown  in  Figures  4.1,  4.2,  the 
only  differences  which  exist  between  them,  are  the 
aerodynamic  mass  parameters  and  nonlinear  derivatives 
changes,  as  referred  to  in  the  (Appendices  A  and  B) . 

Hence,  the  equations  which  represent  the  dynamic 
models  are: 

1.   Pitch  Channel 

^ --q-sfC/v/ w    (gees)  (iv.c.i-i) 

q  -    hi     +  _5HiidOL___ (d€^/sec2)      (lv.C.l-2) 

<i=q  +  K.T1iB-  -|i-  (deglsec)  (iv.c.i-3)  . 

Substituting    the   values    for    aerodynamic    and  mass 
parameter    for   circular    airframe,    the   above   equations 
become : 

71fc&  -  -ZG5292  -Ctf  (IV. G. 1-1) 

q    r-o.om52  p. r+  736.8604 c^  (iv.G.1-2) 
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OC  -    Cj  -  0.9F5M  C^~  O.OiiFH52  p*€  (iv.G.1-3) 
Lateral   Channels 

Vly6Z9^    (CY<.^  +  Cy^-6y)  /W    (g«s)  (IV. C. 2-1) 

«     z    K- viY6  +  a-  p/s*3  -f    UcglseO  (iv.c.2-2) 


r  «~ftr  +  '7*^    (v«+c^.rY*.cSl.W 


(IV. C. 2-3) 

Xxx 

(IV. C. 2-4) 

4m  r    P  Cde^bec)  (iv.c.2-5) 

Substituting    the   values    in    the    above    equations, 
they   yield: 

vwb    -  2.05292  (  CY€.£+  CY<Jy  .£y)  (iv.G.2-1) 

€  r  o.g «5^ (cv£  +  Gyy.<SV)  +0.0^5  p-a-r     (iv.G.2-2) 

V     r-O.01^52  C|.p  +133.333?  (Cnr^  +  Ci^-fy +<^»fe-S<0 

(IV. G. 2-3) 

H.    ANALYSIS  OF  NONLINEAR  CBTT  AUTOPILOTS  FOR  CIRCULAR 

AIRFRAME 

For  purposes  of  analysis,  a  CSMP  program  was  written 
(Appendix  H)  with  the  equations  representing  the 
aerodynamic  models  and  control  laws  of  the  nonlinear  CBTT 
autopilots  for  referred  airframe. 

At  first,  in  order  to  understand  the  effect  of  faster 
response,  in  the  reduction  of  the  effect  of  kinematic  and 
inertial  cross-coupling,  the  pitch  acceleration  error  gain 
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Ka  Figure  4.26  was  reduced  to  K,   =  0.0274  for  slower 
response. 

The  equation,  which  represents  the  acceleration 
compensator  after  the  above  reduction,  is: 

V^-  — " l^C-X)         (IV.H-1) 

5.(-f-  +  0 

Utilizing  method  of  state-space  representation  with 
forcing  term  derivatives,  the  equation  IV.H-1  yields: 

%'   Vh  - 0.95 804  (Hc-X)  (iv.h-2) 

YH  =  -5.V  +q.653-(nic-X)  (iv.h-3) 

Using  the  above  equations,  instead  of  the  equation 
IV. F. 1-2  and  IV. F. 1-3,  a  new  CSMP  program  (Appendix  I)  was 
written. 

Figure  4.28  shows  that  the  achieved  maneuver  plane 
acceleration  has  a  good  response  to  the  first  command  which 
is  applied  at  zero  time.   The  missile  with  the  circular 
airframe  moves  upward  like  a  skid-to-turn  missile,  as  did 
the  elliptical  airframe  missile,  because  the  motion  is  in 
the  desired  maneuver  direction  and  therefore  the  roll 
channel  is  not  commanded. 

The  time  constant  of  the  achieved  maneuver  response  t, 
is  equal  to  0.49  seconds.   The  achieved  maneuver  plane 
response  due  to  the  second  guidance  command  applied  at  3 
seconds,  shown  in  Figure  4.28  is  reacting  differently  to 
the  kinematic  and  inertial  cross-coupling  than  the 
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elliptical  airframe  in  Figure  4.5.   Rather  than  overshoots 
and  undershoots,  a  slowing  transient  starts  at  3.5  seconds. 

Figure  4.29  shows  that  the  overshoot  in  the  body-fixed 
pitch  acceleration,  due  to  the  kinematic  and  inertial 
cross-coupling  during  the  second  command,  is  60.8  percent 
which  is  substantially  more  than  it  was  for  the  elliptical 
airframe  and  occurs  much  sooner  at  3.66  second. 

Figure  4.29  through  4.39  show  the  time  responses  of 

the  nv  ,  a,  8,  q,  p,  r,  4> ,  <Sp,  5  and  5R.   Above  figures 

c 
and  table  III  show  the  characteristic  of  the  responses  of 

the  variables  for  the  circular  airframe,  when  the   K  = 

-0.0274. 

The  referred  results  of  this  study  verify  the  results 

of  [Ref .  2]  . 

I.  EFFECT  OF  INCREASING  PITCH  CHANNEL  SPEED  OF  RESPONSE 
In  order  to  reduce  the  effect  of  kinematic  and 

inertial  cross-coupling  during  the  second  guidance  command, 
the  response  of  the  Pitch  channel  of  the  CBTT  autopilot  for 
the  circular  airframe  was  made  faster  as  shown  in  section 

II.  This  was  accomplished  by  increasing 

the  acceleration  error  gain  Ka  =  -0.0387. 

A2 
Using  the  above  gain  K,   the  equations  represent  the 

A2 
acceleration  compensator  are  the  IV. F. 1-2  and  IV. F. 1-3. 

The  effect  of  the  change  in  KA  on  achieved  body  fixed 
acceleration  (Figure  4.40)  results  in  the  achieved  maneuver 
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Fig.  4.29  Achieved  Body-fixed  Acceleration  (n7  )  vs 

ZB 
Time  (t)  CBTT  of  Circular  airframe;  2  Gees  (0  , 
180°) ,  KA   =  -0.0274 
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Fig.  4.30   Achieved  Body-fixed  Acceleration  (nY  )  vs 

B 
Time  (t);  CBTT  of  Circular  airframe;  2  Gees  (0°, 
180°) ,  Ka  =  -0.0274 
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Fig.  4.31   Angle  of  attack  (a)  vs  Time  (t) 
CBTT  of  Circular  airframe; 
2  Gees  (0°,  180°),  Ka   =  -0.0274 
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Fig.  4.32   Sideslip  angle  U)  vs  Time  (t) 
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Fig.  4.33   Pitch  Body  angular  rate  (q)  vs  Time  (t) 
CBTT  of  Circular  airframe;  2  Gees 
(0  ,  180°) ,  K,  =  -0.0274 
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Fig.  4.34   Yaw  Body  angular  rate  (r)  vs  Time  (t) 
CBTT  of  Circular  airframe; 
2  Gees  (0°,  180°),  K,   =  -0.0274 
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Fig.  4.35   Roll  Body  angular  rate  (P)  vs  Time  (t) 
CBTT  of  Circular  airframes; 
2  Gees  (0°,  180°),  K,   =  -0.0274 
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Fig.  4.36  Roll  angle  U)  vs  Time  (t) 
CBTT  of  Circular  airframe; 
2  Gees  (0°,  180°),  K,  =  -0.0274 
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Fig.  4.37   Pitch  tail  incidence  (5p)  vs  Time  (t) 

CBTT  of  Circular  airframes; 

2  Gees  (0°,  180°),  KA  =  -0.0274 
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Fig.  4.38   Yaw  tail  incidence  (6y)  vs  Time  ( t) 
CBTT  Circular  airframe; 
2  Gees  (0°,  180°),  Ka  =-0.0274 
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Fig.  4.39   Roll  tail  incidence  ( 6R)  vs  Time  (t) 
CBTT  Circular  airframe; 
2  Gees  (0°,  180°) ,  Ka  =-0.0274 
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TABLE  III. 

CBTT  PERFORMANCE  OF  CIRCULAR  AIRFRAME 

VALUE 

TIME(sec) 

CONDITION 

■^  (sec)  ; 
oOershoot% 

x2  (sec) ; 

overshoo 

1.092(gees) 

5.10 

KA 
Al 

0.46  and 

0.52  and 

-3.044 (gees) 

3.060 

% 

% 

1.109(qees) 

4.540 

KA 
A2 

0. 34  and 

0.62  and 

-3.125<gees) 

0.825 

5% 

0.05% 

1.002 (gees) 

5.040 

KA 
A3 

0.34  and 

0.53  and 

-3.125(gees) 

0.825 

5% 

0% 

11.15(deg) 

1.432 

KA 
*1 

0(deg) 

0 

11.489 (deg) 
0.053(deg) 

0.862 

3.656 

KA 
A2 

11.488(deg) 

0.862 

A3 

0(deg) 

0 

4.25(deg) 

3.461 

Kai 

-0.399(deg) 

4.286 

3.633(deg) 

3.458 

A2 

-0.326(deg) 

4.290 

5.114  (deg) 
-0.414(deg) 

3.608 
4.620 

KA 
A3 

21.39 (deg/sec) 

0.285 

KA 

Al 

-8.84(deg/sec) 

3.220 

29.562 (deg/sec) 
-14.323(deg/sec) 

0.261 
3.225 

KA 
A2 

29.562 (deg/sec) 
-17.575(deg/sec) 

0.261 
3.251 

KA 
A3 

278.5(deg/sec) 
-14.166(deg/sec) 

3.394 

4.6 

KA 
Al 

279.487 (deg/sec) 
-14.071 (deg/sec) 

3.4 
4.6 

KA 
A2 

277.574 (deg/sec) 
-13.798 (deg/sec) 

3.4 
4.597 

KA 
A3 

22.66 (deg/sec) 
-1. 456 (deg/sec) 

3.28 
4.56 

KA 
Al 

21.406 (deg/sec) 

3.262 

\ 

-1. 322(deg/sec) 

4.537 

23.159 (deg/sec) 
-1. 364 (deg/sec) 

3.287 
4.818 

\ 
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plane  acceleration  response  (Figure  4.41).   Figure  4.41 
shows  that  the  achieved  maneuver  plane  acceleration 
response  during  the  first  guidance  command  has  improved. 

Table  III  shows  that  increasing  the  KA  resulted  in  a 
decreasing  in  maximum  sideslip  angle  3.66  and  a  slight 
decrease  in  yaw  angular  rates  (maximum  21.4  degrees) . 

J.    EFFECT  OF  INERTIAL  AND  KINEMATIC  CROSS-COUPLING  IN 

PITCH  CHANNEL 

The  cross-couplings,  (-8P)  into  a  and  (rP)  into  q  , 
were  removed  in  order  to  assess  their  effect  on 
performance.   Hence,  the  equations  for  a  and  q  in  the  pitch 
aerodynamic  model  becomes: 

d  -  q  +  k.  Vis      (<Jeg/sec)  (iv.j-d 

9  -5?'3q£dCM  Ijyy  (iv.j-2) 

Replacing  the  values  to  the  parameters  of  the  above 
equations,  they  yield: 

d  -  q  -0.9S54  Ctf  (iv.j-3) 

q    r  ?3S.SG0CH  (iv.j-4) 

where  CN,  C   nonlinear  functions  which  vary  with  a  and 
5   referred  to  in  (Appendix  B) .   The  lateral  aerodynamic 
cross-coupling  was  retained.   For  purpose  of  analysis,  a 
CSMP  program  (Appendix  I)  was  written. 

Comparing  Figures  4.42  and  4.43  with  figures  4.40  and 
4.41  it  is  concluded  that  the  undersirable  transients  are 

gone . 
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Fig.  4.41  Achieved  body-fixed  Acceleration  (n7  )  vs  Time 

ZB 
(t)  CBTT  of  Circular  airframe;  Faster  responding 
Pitch  channel;  2  Gees  (0°,  180°)  K   =  -0.0387 
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Fig.  4.42   Achieved  Inertial  Acceleration  (nz)  vs  Time  (t) 

CBTT  of  Circular  airframe;  Faster  Pitch 
channel;  inertial  and  kinematic  cross-coupling 
into  Pitch  channel  removed;  2  Gees  (0  ,  180  ) , 
KA  =  -0.0387 
A3 
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Table  III  shows  that  the  achieved  maneuver  plane 
acceleration  time  constant  for  the  second  guidance  command 
has  decreased. 

Figures  4.44  through  4.50  show  the  time  responses  for 
a,  B,  q  ,r,  5  ,  &y,    and  <5R  when  the  KA  is  equal  KA=  -0.0387 
and  the  kinematic  and  inertial  cross-coupling  of  the  Pitch 
channel  have  removed. 

K.    CONCLUSIONS 

1.  The  maneuver  plane  responses  have  transients, 
which  may  have  to  be  reduced,  are  caused  by  inertial  and 
kinematic  coupling  between  pitch  and  yaw  dynamics. 
Transients  cause  excessive  overshoots  and  undershoots  in 
achieved  maneuver  plane  acceleration  for  elliptical 
airframe  (Figure  4.6)  and  excessive  slowdown  in  the  speed 
of  response  for  circular  airframe  (Figure  4.28,  4.40). 

2.  The  result  of  the  nonlinear  3-D  performance  study 
verify  the  linear  study  of  chapter  III. 

3.  The  results  of  this  study  verify  the  results  of 
[Ref.  2]. 

Further  analysis  and  conclusion  about  a,  q,  r,  p 
and  6  are  discussed  in  [Ref.  2]. 
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Fig.  4.44   Achieved  Body-fixed  Acceleration  (nv  )  vs  Time 
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(t)  CBTT  of  Circular  airframe;  Faster  pitch 
channel;  inertial  and  kinematic  cross-coupling 
into  pitch  channel  removed;  2  Gees  (0°,  180°) , 


K.  =  -0.0387 
A3 
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Fig.  4.45  Angle  of  Attack  (a)  vs  Time  (t) 

CBTT  of  Circular  airframe;  Faster  Pitch 
channel;  inertial  and  kinematic  cross-coupling 
into  Pitch  channel  removed;  2  Gees 
KA  =  -0.0387 
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Fig.  4.46   Sideslip  angle  ( b)  vs  Time  (t) 

CBTT  of  Circular  airframe;  Faster  Pitch 
channel;  inertial  and  kinematic  cross- 
coupling  into  Pitch  channel  removed; 
2  Gees  (0°,180°),  K.  =  -0.0387 
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Fig.  4.47   Pitch  Body  angular  rate  (q)  vs  Time  (t) 
CBTT  of  Circular  airframe;  Faster  Pitch 
channel;  inertial  and  kinematic  cross-coupling 
into  Pitch  channel  removed;  2  Gees  (0  ,  180  ) , 
K,  =  -0.0387 
A3 
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Fig.  4.48   Yaw  Body  angular  rate  (r)  vs  Time  (t) 

CBTT  of  Circular  airframe;  Faster  Pitch 
channel;  inertial  and  kinematic  cross-coupling 
into  Pitch  channel  removed;  2  Gees  (0°,  180°)  , 
Ka  =  -0.0387 
A3 


184 


e 

o 

/*■ 

e 

IA 

/ 

\ 

r 

o 

\ 

j 

^" 

/ 

go 

/ 

s 

r 

o 

! 

T 

O 

1 

~"\ 

e 

? 

a 

f 

o 
o 

0.0   0.5    1.0    1.9    2.0    2.5    3.0   3.0    4.0    4.5    S.O    5.5    S.O 

TIMElSEC) 

Fig.  4.49   Pitch  tail  incidence  (5  )  vs  Time  (t) 

CBTT  of  Circular  airframe;  Faster  Pitch 
channel) inertial  and  kinematic  cross-coupling 
into  Pitch  channel  removed;  2  Gees  (0  ,  180  )  , 
K,  =  -0.0387 
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Fig.  4.50   Yaw  tail  incidence  (  5y)  vs  Time  (t) 

CBTT  of  Circular  airframe;  Faster  Pitch 
channel;  inertial  and  kinematic  cross-coupling 
into  Pitch  channel  removed;  2  Gees  (0  ,  180  )  , 
Ka  =  -0.0387 
A3 
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V.    MINIMIZATION  OF  THE  KINEMATIC  AND  INERTIAL 
CROSS-COUPLING  EFFECTS 

A.  GENERAL 

The  nonlinear  CBTT  autopilot  for  both  airframes, 
consisted  of  the  uncoupled  autopilots  discussed  in  section 
II,  which  coupled  via  inertial,  kinematic  and  aerodynamic 
couplings . 

The  effects  of  the  cross-coupling  were  discussed  for 
both  airframes  in  section  IV. 

In  this  section  the  kinematic,  inertial  and 
aerodynamic  coupling  effects  will  be  sumarized.   The 
kinematic  and  inertial  effects  in  the  pitch  channel  will  be 
analyzed . 

The  minimization  of  the  above  coupling  and  their 
effects  will  be  studied,  using  feedbacks  of  angle  of  attack 
and  rate  of  angle  of  attacks. 

B.  CROSS-COUPLING  OF  CBTT  AUTOPILOTS 

As  discussed  in  previous  sections,  the  four 
cross-couplings,  between  the  autopilot  channels  for  both 
airframes  are: 

1 .   Aerodynamic  Cross-Coupling 

The  aerodynamic  cross-coupling  is  the  coupling 
between  the  Roll  and  Yaw  nonlinear  channels,  via  the 
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Cn  ,  Co   and  C„   (Figure  4.2)  .   The  major  effect  of  the 

B     6         6 

Y        R 

above  coupling  in  the  elliptical  airframe  response  is  to 
decrease  the  overshoot  in  the  achieved  maneuver  plane 
acceleration  resulting  from  the  second  guidance  command, 
(i.e  2  gee  (180°),  the  other  effects  on  the  missile 
variables  were  small  but  in  a  direction  which  improves 
performance  (i.e  decreased  missile  body  angular  rates,  less 
control  surface  motion,  less  sideslip  variations).   In  the 
circular  airframe,  the  effects  of  the   C^  ,  C   are 

negligible,  the  critical  coupling  coefficient  is  the  C» 

which  effects  the  stability  in  the  coordination  branch, 

which  can  be  improved  by  reducing  the  effect  of  c.   on 

Y 
missile  roll  angular  acceleration.   (equation  IV. C. 2-4) 

2.   Cross-Coupling  Between  Roll-Yaw  Autopilots  for 

Coordination  Motion 

This  cross-coupling  branch,  shown  in  Figure  4.4, 
has  been  added  between  the  yaw  and  roll  control  laws,  in 
order  to  provide  coordinated  motion  between  these  two 
channels.   For  purposes  of  stability  the  coordinating 
branch  gain  Kyp  was  selected  for  the  elliptical  airframe 
Kyp  =  1  and  for  the  circular  airframe  Kyp  =  0.458.   The 
elliptical  airframe  is  more  stable  than  the  circular  with 
the  above  selected   KYP« 
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3.   Kinematic  and  Inertial  Cross-Coupling 

The  kinematic  and  inertial  cross-coupling  exist 
between  the  pitch  and  lateral  nonlinear  aerodynamic  models 
of  the  CBTT  autopilots.   There  exist  the  following 
cross-coupling  paths  between  pitch  and  lateral  channel. 

a.  Kinematic  cross-coupling  path  of  (-Bp)  into  a 
in  the  pitch  aerodynamic  model  (Figure  4.3). 

b.  Kinematic  cross-coupling  path  of  (ap)  into  B 
in  the  lateral  aerodynamic  model  (Figure  4.4). 

c.  Inertial  cross-coupling  path  of  (pr)  into  q  in 
the  lateral  dynamic  model  (Figure  4.4). 

d.  Inertial  cross-coupling  path  of(-qp)  into  r  in 
the   lateral  dynamic  model  (Figure  4.4). 

The  efffects  of  the  inertial  and  kinematic 
coupling  are  the  transients  in  the  achieved  maneuver  plane 
acceleration  response,  and  have  the  form  of  overshoots  and 
undershoots  discussed  in  section  III,  IV.   Transients  may 
also  result  in  a  slower  maneuver  plane  acceleration 
response  and  are  less  pronounced  at  higher  acceleration 
levels. 

Coupling  transients  may  be  reduced  by  increasing 
pitch  stability,  by  decreasing  maximum  roll  rate  (section 
IV) ,  by  improving  the  autopilot  coordination  technique  to 
minimize  sideslip  rate  B  (section  IV) ,  by  increasing  the 
effects  of  stabilizing  lateral  aerodynamic  coupling 
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(negative  C   ) ,  by  increasing  pitch  channel  speed  of 

response  (section  IV)  and  by  increasing  the  pitch  stability 

Cm  with  feedbacks  of  both  angle-of-attack  and  rate 

a 
of  angle-of-attack  for  the  pitch  autopilot. 

C.    KINEMATIC  AND  INERTIAL  CROSS-COUPLING  IN  PITCH 

AUTOPILOT 

The   kinematic   and    inertial    cross-couplings    in    the 
pitch   autopilot,    are   via    the    pitch   angular    rate   q   and    the 
rate   of    angle-of-attack    (a)     (Figure   4.1)). 

The    equations    for    q   and    a    (figure    27)    are: 

9*   ***** fa.*+U.-a)+  ijL     Cv.c.1) 

«  -  °i-   (M/ST-i)  (v.c.2) 

where : 

rP/57.3   =   Inertia  Cross-Coupling 

-6P/57.3  =   Kinematic  Cross-Coupling 

Removing  C   from  the  parenthesis  in  (V.C.I)  and 
a 
combining  with  (V.C.2)  to  eliminate  q,  the  (V.C.I)  results 

in: 


T       U>ia    '   5T.3   >>w3  ' 


(V.C.3) 
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Let: 

Iyy 

and 

u2=  -— £- 

Taking  the  Laplace  transform  of  (V.C.3)  and  solving 
for  a,  the  equation  becomes: 

-j        f(rP/5T.3)-S(M)        X  \   , 
OC= r •   1^J- fK2-d(>  )    (V.C.4) 

Increasing  the  magnitude  of  Cm  will  decrease  both  K^ 

a 
and  1/K,  resulting  in  the  reduction  of  cross-coupling 

effects  on  the  angle-of-attack  (a)  and  the  influence  of  5_ 

in  the  a . 

In  the  section  III  the  equation  of  n_  was  derived: 

B 

Mj6r  (-qrf/w)  •(<«,•$>+ Gb*«)  (v-c-5> 

which  involve  (a).   Therefore,  the  reduction  of  inertial 

and  kinematic  effects  on  (a) ,  results  in  the  reduction  of 

the  effects  in  the  achieved  body-fixed  acceleration  (n7  ) . 

ZB 
In  the  nonlinear  studies  the  equation  for  the  achieved 

inertial  acceleration  was  derived  as: 

^l'   %  'C05<J>+71y6.5lT)<b  (V.C.6) 

It  is  apparent  from  the  equations  (V.C.4) ,  (V.C.5)  and 

(V.C.6) ,  that  an  increasing  of  C   results  in  the  reduction 

a 
of  kinematic  and  inertial  effects  (i.e  transients  in  the 

response)  in  the  achieved  inertial  acceleration  (n~) ,  which 
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it  is  desired  in  order  to  meet  the  requirement  of  time 
constants  (?,,  t^)  and  percents  of  overshoots  and 
undershoot,  as  referred  in  Appendix  C. 

D.    MINIMIZATION  OF  THE  KINEMATIC  AND  INERTIAL 

CROSS-COUPLING  EFFECTS  IN  ELLIPTICAL  AIRFRAME 
As  discussed  in  the  nonlinear  studies,  using  the  pitch 
angular  rate  q  as  feedback  and  under  the  influence  of  the 
inertial  and  kinematic  effects,  the  achieved  inertial 
acceleration  (nz)  for  both  airframes  doesn't  meet  the 
requirements  of  the  time  constants,  overshoot  and 
undershoot,  denoted  for  the  linear  autopilots,  referred  to 
in  Appendix  C. 

Figure  4.5  and  Table  II  show  that  the  achieved 
inertial  acceleration  (n„)  for  elliptical  airframe  with 
feedback  the  pitch  angular  rate  (q)  has  overshoot  11.25% 
and  t1  =  0.48  sec,  *2  =  0-4l  sec» 

Figures  4.40  and  Table  III  illustrate  that  the 
achieved  inertial  acceleration  (n„)  for  circular  airframe, 
with  feedback  the  (q)   has  overshoot  5%  and  t=-.46  sec, 
x     =  -.62  sec.   Hence  the  (nz)  for  both  airframes  is  not 
acceptable. 

In  order  to  minimize  the  inertial  and  kinematic 
cross-coupling  effect,  namely  to  have  smoother  transient 
response  of  n„  which  will  meet  the  requirements  of  time 
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constants  and  overshoot,  the  magnitude  of  C   was  increased, 

a 
utilizing  as  feedback  the  angle-of-attack  and  its 

derivative. 

In  order  to  increase  the  magnitude  of  C„   the 

3  m 

a 

stability  augmentation  theory  was  used. 

Considered  as  stability  derivative  augmentation  a— *6 

P 

and  a-*(S   increase  the  magnitude  of  Ma  and  M^  directly,  as 

referred  in  [Ref.  6],  the  increases  in  the  Ma   is 

essentially  an  increase  in  C    and  hence  the  inertial  and 

a 
kinematic  effects  on  the  angle-of-attack  are  reduced. 

In  the  pitch  control  law  (Figure  3.1)  a  SAS  exists 

using  as  feedback  the  pitch  angular  rate  (q) .   In  order  to 

design  the  SAS  with  the  angle  of  attack  (a)  as  feedback, 

the  root  locus  and  lead-lag  compensation  techniques  were 

used . 

1 .   Stability  Augmentation  System  (SAS)  With  Feedback 

of  the  Pitch  Angular  Rate  (g) . 

In  Figure  5.1,  the  stability  augmentation  system 
is  shown,  with  feedback  of  pitch  angular  rate  q  where  the 
transfer  function  of  the  q/<5   for  the  uncoupled  pitch 
dynamic  model,  given  in  section  II. C.   The  closed-loop 
transfer  function  of  the  above  SAS  is: 

G(S)z _ (V.D.l-l) 

1.,  k(-H)(5*3)(-*H.35-5-*.?H*) 

S(5?.3)(^  +  l)(s240.1SJS-li.2T9) 
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Manipulating    equation   V.D.1-1,    it    results: 
K  (-236.183?  S2-3W3fe3-3  -MHl-8862) 


fa)- 


0.0053- SM+1.00O9-5Vo.af-i1-lJ.28S+K(-5.l69S-^.32U-7.?H8) 

The  poles  and  zeros  of  the  SAS  are  as  follows: 

a.  P±   =  -188.43 

b.  P2  =  -3.45 

c.  P3  =  3.269 

d.  P4  =  0.0 
and 

e.  Z,  =  -8.0 

f.  Z2  =  -0.186 

The  other  two  zeros  are  at  infinity. 
The  root  locus  of  the  system  is  shown  in  Figure 
(5.2). 

The  K  was  selecting  at  the  value  of  -3.07,  because 
at  this  value,  the  system  has  two  dominant  roots  at  P,,2  = 
-9.212 ±j7. 606  which  dominate  the  response  of  the  system, 
having  the  following  characteristics: 

a.  8:W}-Wi|--tf8 


where : 


d.  LUt,  =  A  =  K.203  (voui/s*c) 
d 


damping  constant 
frequency  of  system 


(u 


0) 


natural  undamped  frequency 
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RO0LOC/BIPOP  ANALYSIS  OF  CBTT  AUTOPILOT 
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Fig.  5.2   Root  locus  of  SAS  with  q/6D/  of  pitch 
channel;  CBTT  of  Elliptical  airframe 
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2.   SAS  with  Feedback  of  the  Angle  of  Attack  (a)  and 

rate  of  angle  of  attack  (a) 

The  SAS  with  feedback  the  pitch  angular  rate  q 
(Figure  5.1)  and  K  =  -3.07,  is  stable  and  the  dominant 
roots  gives  damping  ratio  5  =  0.73  in  the  response  of  the 
system. 

To  minimize  the  coupling  effects  the  angle  of 
attack  (a)  was  selected  as  feedback  in  the  previous  SAS. 
It  is  desired  the  SAS  to  have  the  same  behavior  as  before, 
for  purposes  of  stability  and  response  of  whole  the  pitch 
channel . 

Pole  zero  cancellation  technique  is  used,  to 
obtain  the  same  poles  and  zeros,  as  in  the  case  with 
feedback  the  pitch  angular  rate  q. 

For  above  purposes,  a  compensator  was  used  in  the 
SAS,  as  shown  in  Figure  5.3. 

The  transfer  function  of  the  compensator  is: 

GW-  -O.02.S-HI.35*  (V.D.2-1) 

The  closed  loop  transfer  function  of  the  SAS  is  as 
follows: 

\-0.0*S-m.35?  )    {        5    ly^-gj  V324  0tJglS-iJ.2?9  ) 


G<CS)  = 


il(i   \/-MUS?-5-KFHg\    MjlOU   5*»*     U-O-02-S-*H.35T       \ 


(V.D.2-2) 


where  the  transfer  function   a  _   -0-Q3-S  -m.357  ±s 


derived  in  the  section  II. C. 3. 
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Rearranging    and   manipulating    the    equation 
(V.D.2-2),    it    becomes: 

K  (5.9Z315  +  1229.381?+  100225. MM9-S  4-  1$2J5.0<j) 


tor 


and 


-0.0001-S5+  (X2215  SS-<H.399  55-<U?9  S2-  466.q66S  +  K  (Q.103M  S3+2IM.6  s'+WOU  S+3IM) 

(V.D.2-3) 
The  poles  and  zeros  of  the  system  are: 
P-L  =  -1877.931 
P2  =  -207.391 
P3  =  -3.447 
P4  =  -3.271 
P5  =  0.0 

Zx  =  -2067.838 

Z2  =  -8.0 

Z3  =  -0.186 

The  other  two  zeros  are  at  infinity. 

The  root  locus  of  the  system  is  shown  in  figure 


5.4. 


Selecting  the  K  =  -3.07,  the  system  has  the  same 
dominant  roots  P-j^  2  =  -8.249  ±j7 .  577,  as  in  section  V.D.I. 
3.   Acceleration  Loop 

Considering  the  achieved  acceleration  loop  in  the 
Pitch  channel,  the  modification  of  the  SAS,  effects  the 
transmission  gain  of  the  above  loop. 
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In  order  to  compensate  the  effect  in  the 
transmission  path,  a  compensator  was  added  as  shown  in 
Figure  5.5. 

The  transfer  function  of  the  compensator  is  the 

inverse  of  the  V.D.2-1: 

r   -  N   1    -0.02.S-iH.35? 

Modifying  in  this  way,  the  transmission  path  of 
the  acceleration  loop,  the  gain  of  path,  is  the  same  as  in 
the  case  of  par.  V.D.I. 

4.   Measuring  of  the  Minimization  in  the  Kinematic  and 

Inertial  Cross-Coupling  Effect 

In  order  to  measure  the  minimization  a  CSMP 
program  (Appendix  J)  was  written  using  the  equations  of  the 
nonlinear  aerodynamic  models,  the  equations  of  nonlinear 
lateral  control  law  and  the  equations  of  the  modified  pitch 
control  law. 

The  equations  of  the  modified  pitch  control  law 
are  the  following: 

a.  Acceleration  lag  network 

X:  -150X4150  H*b  (II. C. 4-2) 

b.  Acceleration  second  order  compensator  (Figure 

5*5)  V         (-0.008).  (-0.03- S-MJ.35?) 

(U1VX)  "   ($+fl(-U35*5-?w«0    (V-D-4"1) 
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Utilizing    inverse   Laplace    transformation   and 
rearranging    the    above    equation,    one   has: 

S+U8G?Y+;ms  V  =  -0.000232 (Lffilc-X)-a4M-(4-ilVX)  (V.D.4-2) 
Using    space    representation    of    a    system,    in   which 
the    forcing    function    involves   derivative 
[Ref.    5,    pp   675-678],    one    obtains: 

Y/f  Z  X2-0.OOO25?5-T1c  +0.000232  -X  (V.D.4-3) 

W2  z  -1.1188?  W,  -  6.1S6?-V^  -0.53i2Hc+0.^m-X         (V.  D .  4-4 ) 

V  =  WA  (V.D.4-5) 


(V.D.4-5a) 


c.  Compensator    of   Actuator 

ftc         .     K(-r44>>nf.357-5-?.?118) 
(V-a|s*3)  "  S  (-0.02S-M1.35?) 

Using    the   method    for    state    representation    and 

rearranging    the    equation  V.D.4-5   yields: 

<Jp  Z(5p2+K-  532332 Xl03-(M-al5T.3)  (V.D.4-6) 

S?zz  -20^ZSS^-  K.  1100.SS6X10C(N-al57-3)  (V.D.4-7) 

S^z  <5>f  -  K-  2Sg.48(Y-a|57.3)  (V.D.4-8) 

d.  Equation    of   Actuator 

Jp, -18S-M- Jp4  -10*95.3  -opc  (IV. B.  1-4) 

Selecting  K=3.07,  for  the  compensator  of  SAS,  with 
feedback  (a)  the  time  constant  t..  (i.e  time  constant  of 
achieved  maneuver  acceleration  due  to  the  first  command)  is 
0.45  sec,  the  time  constant  x   (for  second  command)  is  0.41 
sec  and  the  undershoot  is  20.73  percent  as  shown  in  figure 
5.6. 
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The  above  response  of  the  inertial  achieved 
maneuver  acceleration  does  not  have,  smooth  transient, 
because  of  the  large  undershoot. 

Due  to  the  above  reason,  another  K  was  selected 
equal  to  -19.37,  for  minimization  of  the  coupling  effects. 

With  this  selected  K  the  system  has  the  following 
poles,  shown  in  figure  5.4 

Px,2  =  -98±j96.239 

P3   =  -1879.06 

P4    =  -8.596 

P5    =  -0.189 
With  the  same  zeros  as  section  V.D.2,  the 
dominant  poles  of  the  system  are  P-^2  and  tne 
characteristics  of  the  response  are: 

a.  6  =    44.24 

b.  c   =    0.716 

c.  ui  =    93.775    rad/sec 

d.  u  =   134.412    rad/sec 

It    is    shown    in    known    (Ref.    6)    that    the    ">n   can   be 
expressed    as: 


(V.D.4-9) 
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and 


&   m,2  a  (V.D.4-10) 


Ma=  C„,a.-f-.p.M'.A 


where: 

y:   specific  weight  of  air  1.4 

I:   moment  of  inertia 

P:   Local  pressure 

M:   Mach  number 

A:   Area 

Comparing  the  u  =  134.412  for  the  K  =  -19.37,  it 

is  apparent  that  it  is  much  larger  than  the  o>n  =  11.203  for 

K  =  -3.07  and  hence  according  to  the  equations  V.D.4-9  and 

V.D.4-10,  the  Cm  increases  essentially,  resulting  in 

a 
minimization  of  kinematic  and  inertial  cross  coupling 

effects . 

Figure  5.7  shows,  that  with  K  =  -19.37  and 
feedback  in  SAS  the  (a),  the  achieved  inertial  maneuver 
acceleration  has  t  =  -.49  sec,  t2  =  -.38  sec  and  overshoot 
about  0%. 

The  above  analysis  illustrates  that  the  selection 
of  K  =  -19.37,  instead  of  K  =  -3.07,  keeps  the  damping 
ratio  about  constant  c  =  0.71  (c  =  0.73  for  K  =-3.07)  and 

increases  the  w  .  which  essentially  increases  the  C   , 

n'  *  m 

a 

resulting  in  very  smooth  response  of  the  achieved  maneuver 
acceleration. 
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Using  the  same  particular  compensators  in  the 
pitch  channel  shown  in  Figure  5.5  and  feedback  the  rate  of 
angle  of  attack  (a),  the  Pitch  channel  becomes 
uncontrolable  for  both  selection  of  K  =  -3.07  and  K  = 
-19.37,  because  the  pitch  tail  incidence  exceeds  the  limits 
of  ±10°. 

Also,  selecting  the  K  =  -19.37  in  SAS  of  Figure 
5.1  with  feedback  the  pitch  angular  rate  q  the  Pitch 
channel  becomes  uncontrolable  again,  due  to  the  exceeding 
of  the  limits  in  5p. 

Therefore,  for  minimization  of  the  kinematic  and 

inertial  effects  (i.e  improving  of  the  transient  response 

of  the  achieved  maneuver  acceleration)  for  this  particular 

airframe  and  flight  conditions,  the  Cm  was  increased, 

a 
using  the  aforementioned  compensators  in  the  pitch  control 

law  (Figure  5.5)  and  feedback  the  angle-of-attack  (a). 

Comparing  the  achieved  inertial  acceleration  (n~) 

with  feedback  (q)  in  the  Pitch  channel  (Figure  4.41)  and 

the  achieved  inertial  acceleration  (nz)  with  feedback  (a) 

and  K=  -19.37  (Figure  5.7),  it  is  concluded  that  for  these 

particular  airframes  and  flight  conditions  the  feedback  of 

angle-of-attack  (a)  is  more  desirable,  than  the  feedback  of 

(q)  and  (a) ,  because  in  the  achieved  inertial  acceleration 

has  been  reduced  the  kinematic  and  inertial  effects, 

meeting  with  this  way  the  requirement  of  the  time  constants 

and  overshoot,  denoted  in  the  linear  studies. 
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E.    MINIMIZATION  OF  KINEMATIC  AND  INERTIAL  COUPLING  IN 

ELLIPTICAL  AIRFRAME 

Making  a  comparison  in  Figure  5.8,  which  shows  the 
kinematic  cross-coupling  (-pB)  into  a,  in  the  pitch  channel 
with  feedback  pitch  angular  rate  (q)  and  Figure  5.10  which 
shows  the  -p&  in  the  modified  pitch  channel  (Figure  5-5) 
with  feedback  the  angle-of-attack,  it  is  concluded  (Table 
IV)  that  the  kinematic  coupling  is  increased  from 
(-0°. 988-5°. 877)  to  (-2°. 124-7°. 383) ,  when  the  kinematic 
coupling  effects  are  reduced  by  the  modified  pitch  channel. 

Figure  5.9  shows  the  inertial  coupling  (pr)  into  q  in 
the  nonlinear  Pitch  channel,  with  feedback  the  pitch 
angular  rate  (q)  . 

Figure  5.11  shows  the  inertial  coupling  for  the 
modified  Pitch  autopilot  with  feedback  the  angle  of  attack 
(a)  (Figure  5.5)  . 

Comparing  the  above  figures,  it  is  concluded  (Table 
IV)  that  the  inertial  coupling  is  increased  from 
(-1°. 692-168°. 317)  to  (-20°. 373  -  181°. 415),  when  the 
inertial  coupling  effects  are  reduced  by  the  modified  pitch 
channel  discussed  in  previous  section. 

Hence,  the  minimization  of  the  kinematic  and  inertial 
effect  does  not  mean  necessarily  minimization  of  the 
cross-couplings  themselves. 
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For  minimization  of  the  kinematic  and  inertial  cross 
coupling,  the  same  modified  nonlinear  pitch  control  law  was 
used  as  that  of  section  IV. D  (Figure  5.5),  except  the 
acceleration  compensator  which  was  removed  from  the 
transmission  path  of  the  nQ  as  shown  in  Figure  5.12. 

For  purposes  of  analysis  a  CSMP  program  (Appendix  J) 
was  written,  using  the  following  equation  for  the  nonlinear 
pitch  control  law  and  the  same  equation  for  lateral  control 
law  and  dynamic  models,  as  in  section  IV. B  and  IV. C. 

X    r   -  150X  +  i50-T)*8  (iv. b. l-l) 

S^    -   -6-V  -0.5305  7\c  +  0.48- X  (iv.b.i-2) 

Sty  -  Sp2  +  K.  532.382X403-  (Y-a|5?.$)  (v.d.4-6) 

<Sp2  -  -2O^.95.(S"p2-K.HO0.m^06(V-al57.3)     (v.D.4-7) 
£Pc=  5?^  -K-  25S-MS  (Y-a|5?.3)  (v.D.4-8) 

Sp  -  -18*.M-5p  +  40}953-£pc  (iv. b. 1-4) 

where  K  =  19.37. 

Figure  5.13  shows  the  minimized  kinematic  coupling  and 
Figure  5.14  shows  the  minimized  inertial  coupling.   Figures 
5.13,  5.14  and  table  IV  illustrate  that  the  minimization  of 
the  kinematic  and  inertial  coupling  is  264%  and  195.26% 
respectively  from  its  values  in  the  nonlinear  CBTT 
autopilot  for  elliptical  airframe.   Figure  5.15  shows  the 
achieved  inertial  acceleration  (nz) ,  with  the  Pitch  channel 
modified  for  minimization  of  kinematic  and  inertial 
coupling  (Figure  5.12).   From  this  figure  it  is  concluded 
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TABLE  IV 


KINEMATIC  AND  INERTIAL  CROSS-COUPLING  FOR 
ELLIPTICAL  AIRFRAME 

CONDITION  INERTIAL  (deg/sec2)   KINEMATIC  (deg/sec2) 


MIN        MAX         MIN         MAX 
Elliptical  airframe 
with  cross- 
coupling  -1.692     168.847     -0.9881      5.877 
(section  IV. C) 


Elliptical  ideal 

airframe  (section    -1.698     167.37     -2.124      7.383 

IV. E) 


Elliptical  airframe 

with  modified  Pitch 

autopilot  -20.373    181.415    -1.88       8.644 

K  =  -19.37 

(section  V.l) 


Modified  Pitch 

autopilot  K  =  -30.4   -7.263     169.371     -1.956       8.711 

(section  V.D. ) 


Modified  Pitch 

autopilot  for        -0.0577    57.5322    -0.523      1.363 

minimization  of 

cross  couplings 

(section  V.E) 

K  =  -19.37 


Note:  K  is  the  gain  of  the  added  compensator  in  SAS  of  the 
Pitch  autopilot. 
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that  the  kinematic  and  inertial  effects  (i.e  transient  in 
achieved  acceleration  response)  is  increased. 

F.    MINIMIZATION  OF  THE  KINEMATIC  AND  INERTIAL 
CROSS-COUPLING  EFFECTS  IN  CIRCULAR  AIRFRAME 
As  shown  in  Figures  4.28,  4.40  and  4.42  the  achieved 
maneuver  acceleration  in  the  circular  airframe  has  very 
smooth  response,  in  regards  to  the  overshoot  or  undershoot 
on  it. 

The  effect  of  the  inertial  and  kinematic 
cross-coupling  in  the  nonlinear  aerodynamic  model  is  that 
the  response  of  the  inertial  achieved  acceleration  has  a 
slowing  transient  starts  at  3.5  seconds,  as  discussed  in 
section  IV. G.   Hence,  minimization  of  the  effects  means  in 
the  circular  airframe  improvement  of  time  constant  during 
the  second  guidance  command. 

1.   Stability  Augmentation  System  (SAS)  with  Feedback 
the  q. 

Using  the  same  method  for  minimization,  as  in 
section  V.D,  the  SAS  of  the  pitch  channel  for  the  circular 
airframe  with  feedback  q  is,  as  shown  in  Figure  5.16.   The 
transfer  function  of  the  q/5D  is  given  in  section  II. C. 

The  closed-loop  transfer  function  of  the  above  SAS 
is: 
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(lis*  A)  5?-3   (  "  5M°»W  -  6.  W5> 
2W  = — k  affy  Iiii**3 — LI (v.F.i-i) 

K-  (  csS+1)  (ST-3)  (-53JQS8S  ^.qH5) 

Rearranging  equation  V.F.1-1,  one  obtains: 
K  {20.tt'?±30f00.2SS2+Q038n6S+23tJ2.?)         

(V.F.1-2) 
The  poles  and  zeros  of  the  SAS  are:  (1) 
P,=  -188.438;  (2)  P2  3  =  -0.0739^6.657;  (3)  P4  =  -0.143; 
(4)  Z-,  =  -6.55  and  (5)  Z2  =  -0.118.   The  other  two  zeros 
are  at  infinity.   The  root  locus  of  the  SAS  is  shown  in 
Figure  5.17.   The  selected  value  of  K  is  -15.6  and  the 
dominant  roots  are  P,  2  =  -10. 736±j7 .  536  (Figure  5.17). 
With  those  dominant  roots,  the  response  of  the  system  has  9 
equals  to  35.06  degrees;  C  equals  to  0.818;  w  equals  to 
7.545  rad/sec  and  <*>   equals  to  13.125  rad/sec. 
2.   SAS  with  Feedback  of  the  a 

For  minimization  of  inertial  and  kinematic  effect 
in  the  circular  airframe,  the  same  technique  was  used  (i.e 
pole-zero  cancellation  technique)  as  in  section  V.D.2  and 
the  resulted  pitch  control  law  with  feedback  of 
angle-of-attack  (a),  is  shown  in  Figure  5.18. 

The  transfer  function  of  the  added  compensator  in 
the  SAS  is  given  by  equation  V.D.2-2. 

222  (V.D.2-2) 


REAL  RX1S [UNITS  PER  INCH)  =  20.0000 

IMflG   AXIS  [UNITS   PER    1NCHI     =    10.0000  + 

R00L0C/BIP0P  ANALYSIS  OF  CBTT  AUTOPILOT 
CIRCULAR  AIRFRAME 


Fig.  5.17   Root  locus  of  SAS  with  q/ 6  ;  CBTT  of 
Circular  airframe 
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The   closed    loop    transfer    function   of    the    referred 
before    system    is    given    by   equation  V.D.2-3: 

/K  fas?*1 )  \   As9.i-S-6.93  \    /   5T3     \     /-0.039-S-59.1        \ 
V  (<Hm3 +^    J    V-0.O39S-59.1  J   [j$tt  +  \    )     \  S%0>ms+W.lf~  I 

K.(^+h)5?.3-(-59.1S-6.99)  (-0.039 5- 5S.1) 

(V.D.2-3) 
where  the  transfer  function  &/ &n    is  given  in  the 
section  II. C. 3,  rearranging  and  minor  manipulating  the 
above  equation,  one  obtains  equation  V.D.2-4: 


60= 


H  (20.38- &  303QQ.2S  S*+  Z0S821. 6S-h  236*2, 3  ) 


~0.OOlS*-2MlSH-4m.3MS  -229.7  5-f*3<f9-5-i630.4+fc  ((U6  5V5a5.«S+355»+^i3) 

(V.D.2-4) 

The  poles  and  zeros  of  the  modified  SAS  are, 

P1   =  -1499.623,  P2  =  -188.53,  P3,4  =  -0 . 0739 ±j6 . 656 ,  P5  = 

-0.142,  Z1  =  -1499.6,  Z2  =  -6.55,  Z3  =  -0.1182  and  all 

other  zeros  at  infinity. 

The  root  locus  of  the  system  is  shown  in  figure 

5.19,  it  is  apparent  from  above  Figure,  that  for  k=-15.6 

the  dominant  roots  are  the  same  as  those  in  section  V.F.I. 

Considering  the  transmission  of  acceleration  n 

c 

path,  as  in  the  elliptical  airframe,  a  compensator  was 
added  in  the  path,  for  compensation  of  the  effect  in  it, 
which  is  created  by  the  added  compensator  in  the  SAS,  as 
shown  in  Figure  5.20. 
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REAL  RXIS  UNITS  PER    INCH!    =  20.0000 
IHRG  flXISCUNITS  PER    INCHI    =    10.0000 

RaOLOC/BIPOP  ANALYSIS  OF 
CIRCULAR  AIRFRAME 


CBTT  AUTOPILOT 


Fig.  5.19   Root  locus  of  SAS  with  a/5 
CBTT  of  Circular  airframe 
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3.   Measuring  of  the  Minimization  in  the  Kinematic  and 
Inertial  Cross-Coupling  Effect 
A  CSMP  program  (Appendix  K)  was  written  using 
equations  of  nonlinear  lateral  control  law  (section 
IV. G. 2),  of  the  aerodynamic  models  (section  IV. 4)  and  the 
following  equations  of  the  modified  nonlinear  pitch  control 
law: 

a.  Acceleration  lag  Compensator 

X  r  -150X  +150  Yfe6  (II. C. 4) 

b.  Acceleration  Second-Order  Compensator  (Figure 

3-20)  . 

_^_  r  K.  U^-^J  rOggjg^l)(v,r.3-i) 
i\-X)         5  (|-+l).  (-59.10S-S- 6.38^5) 

Utilizing  inverse  Laplace  transformation  and 
rearranging  the  above  equation,  yields 

Y3  +5.1183  Y  +  0.591  Y  r-k. 0.022*8  (\-  X  )- 3U.9?(vX)-  S^X) 

(V.F.3-2) 
Using  state-space  representation  referred  to  in 
(Ref.  5:  pp  675-678),  one  obtains: 

H  -  W2-  KA-  0.022 «(VX)  (v.f.3-3) 

Vf2-  \M3-  kA-  3M.S50  (t^-X^  (v.f.3-4) 

%z  -0.5eU.W2-5.118.^34.\(v1?3.3!)(vx)  (v. p. 3-5) 
V3  =   YV<l  (v.F.3-6) 
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c.  Compensator    of   Actuator    in    SAS 

j>c  _    K.(^5-H)(-59J0ggS-6939)    (y.F.3-7) 

(V  Ot/5f3)         (caW+4)"  (-0.03945- 59.1(H) 

Utilizing,  the  aforementioned  method  for  state 
representation  and  rearranging  the  equation  V.F.3-7,  it 
yields : 

$?{  -  <%  4  K  (w<36.0*33).  (Y3- CCJ57.3)   ( V .  F .  3-8 ) 
$?2  --mWS-fy  -  1500  <5p2-  V.^Xtofyatett^V.  F.  3-9  ) 

<Tpc=  Jft  -  K.  (52.W)  (Y3-a/573)     (V.F.  3-10) 

d.  Equation  of  Autuator 

(IV. G. 3-8) 


£p  =  -i*M  -5p+  *0795.3.Sf> 


Selecting  the  K  =  -15.6  and  K*  =  -0.0387 

A2 
(Figure  4.26)  and  feedback  the  pitch  angular  rate  (q) ,  the 

achieved  inertial  acceleration  n„ ,  has  T-,  =  -.34  sec,  T2  = 
0.62  sec  and  overshoot  5%  (Figure  4.40,  Table  III). 

This  response  of  the  n„  is  unacceptable 
because  it  doesn't  meet  the  requirement  on  the  time 
constant,  referred  to  in  Appendix  C. 

Selecting  K  =  -15.6  and  KA  =  -0.0387  (figure 
5.20)  and  feedback  the  angle  of  attack  (a)  in  the  pitch 
channel  the  achieved  inertial  acceleration  nz  (Figure  5.21) 
has  t  =  -.35  sec,  t  =  0.6  sec  and  overshoot  7%,  which  is 
again  unacceptable. 
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It    is    known    (Ref.    6)    that: 


Xx*  |f>M2-A 


(V.G.3-9) 


Hence,  in  order  to  reduce  the  t   it  is  needed 

to  increase  the  Cm  .   Selecting  K  =  -70.88  and  KA  =  -0.0387 

a 
(Figure  5.20)  with  feedback  the  angle  of  attack  (a),  the 

dominant  roots  of  the  SAS  (Figure  5.18)  are  P-,,2  = 

-90. 55±j87.95  and  the  characteristics  of  the  system 

response  are  9  =  44.16  degrees,  c  =  0.79  ,  a  =  87.932 

rad/sec  and  u  =  126.212  rad/sec.   With  that  selected  k  the 

z  decreases  from  0.818  (when  K=  -15.6)  to  0.72  and  the  « 

increases  from  13.125  rad/sec  (when  K  =  -15.6)  to  126.212 

rad/sec. 

Figure  5.22  shows,  the  achieved  inertial 

acceleration  n7 ,  when  K  =  -70.88,  Ka  =  -0.0387  and  feedback 

L  A2 

the  (a).   The  t.  =  0.33  sec,  t  =  0.48  sec  and  overshoot 

14.33  %.   It  is  apparent  that  the  t   decreases,  which  is 

desirable  but  the  overshoot  increases,  because  the  damping 

ratio  x.   decreases. 

In  order  to  reduce  the  overshoot  and  to  keep 
the  time  constants  below  the  0.5  sec,  the  KA  reduces  to 
-0.0274  (Figure  5.20)  and  the  gain  Kyp  (Figure  4.27)  of  the 
coordination  branch  in  the  nonlinear  lateral  control  law 
increases  from  0.458  to  1.0. 

Using  Kvp  =  1.0  (Figure  4.27)  in  the 
coordination  branch  in  the  lateral  control  law,  KA 
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LEGENO 
NZ  □ 


u.eo         s.«o         e.ao 


'0.00  O.SO  1.20  1.80  Z'.U  9.00       x      3'. 

TIMEQSeC) 

Fig.  5.25   Achieved  Inertial  Acceleration  (n  )  vs  Time  (t) 


CBTT  of  Circular  airframe;  K  =  -70.88, 
KA  =  -0.0274,  Kvn  =  1.0 
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=  -0.0274  in  the  acceleration  compensator  and  K  =  -70.88  in 
the  added  compensator  in  the  SAS  with  feedback  the 
angle-of-attack  (a)  (Figure  5.20),  the  obtained  achieved 
maneuver  acceleration  has  T,  =  -.43  sec,  t   =  -.37  sec 
overshoot  8.9%,  as  shown  in  figure  5.25  and  table  V. 

Using  in  the  pitch  channel  with  the 
aforementioned  particular  compensators,  as  feedbback  the 
rate  of  angle  of  attack  (a)  and  making  all  the  above 
selection  of  KA,  K  and  Kvp,  the  whole  system  becomes 
uncontrolable  because  the  required  tail  incidence  angle 
exceeds  the  limits  of  ±10°. 

Also  using  the  selections  of  K,  K,,  Kyp  and 
feedback  pitch  angular  rate  q  (Figures  4.26,  4.27),  again 
the  whole  system  becomes  uncontrolable,  because  the 
required  tail  incidence  angle  exceeds  the  limit  of  ±10  . 

The  above  analysis  illustrates  that  the 
angle-of-attack  is  more  desirable  feedback  than  the  (q)  and 
(a) ,  in  the  nonlinear  pitch  control  law  for  these 
particular  circular  airframe  and  flight  conditions,  because 
it  reduces  the  kinematic  and  inertial  coupling  effect,  such 
that  the  achieved  inertial  acceleration  (nz)  meets  the 
requirements  of  the  time  constants  and  overshoot,  referred 
to  in  the  linear  studies. 
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TABLE  V 


CHARACTERISTIC  OF  THE  INERTIAL  ACHIEVED  ACCELERATION  (nz) 
FOR  CIRCULAR  AIRFRAME 


CONDITION 


t,  (sec)     t 2  (sec)    overshoot 


0.6 

7 

0.48 

14.33 

0.51 

8.6 

KA  =  0.0387,  K  =  -15.6  0.35 

KA  =  0.0387,  K  =  -70.88  0.33 

KA  =  0.0274,  K  =  -15.6  0.435 

KA  =  0.0274,  K  =  -70.88 

Kyp=  1  0.43  0.37        8.9 

Note:   K.:  pitch  autopilot  acceleration  gain  in  acceleration 
compensator 
t,:  63  percent  time  constant  of  achieved  maneuver 

plane  acceleration  due  to  first  command 
?2:  63  percent  time  constant  of  achieved  maneuver 
plane  acceleration  due  to  second  command 
Kyp:  gain  of  the  coordination  branch 

K:  gain  of  the  added  compensator  is  the  SAS  of  pitch 
autopilot 
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G.    MINIMIZATION  OF  KINEMATIC  AND  INERTIAL  COUPLING  IN 

CIRCULAR  AIRFRAME 

As  in  the  elliptical  airframe  comparing  the  values  of 
the  inertial  and  kinematic  cross-coupling  for  the  circular 
airframe,  namely,  the  (P.r)  and  (-p$) ,  referred  to  in  table 
VI,  it  is  concluded  for  this  particular  analysis  that  as 
the  coupling  effects  in  the  nonlinear  Pitch  autopilot  are 
reduced,  the  cross-couplings  themselves  are  increased. 

In  order  to  minimize  the  cross-couplings,  the  pitch 
control  law  was  modified,  as  shown  in  Figure  5.26. 

For  purposes  of  analysis,  a  CSMP  program  (Appendix  K) 
was  written,  using  the  following  equations  for  the 
nonlinear  modified  pitch  control  law,  lateral  control  law 
and  dynamic  models,  as  in  sections  IV. G  and  IV. H. 

1-  X  r  -150X  +150  -^B  (ii. c. 4) 

2.  W^r  W2-V0.022S8(V>^  (v.f.3-3) 

3.  W2r  W3-  Ka'34.850(t|c-X)  (v.f.3-4) 
4-  ^3--0.59lW2-5.^8TW3+KAH?3385(vx)  (v.f.4-5) 

5.  V3=  Y/i  (V.F.4-6) 

6.  Jpc=  -a4M3-Spc-0.3^0S8(Va/5?.3)-2230^-(%-a/5>.3) 

(V.F.4-7) 
which  is  the  equation  IV. G. 1-4,  except  of  the 
feedback,  which  is  the  angle  of  attack  (a)  instead  of  the 
pitch  angular  rate  (q)  . 

7.  <$j>  r  H8$.M-£p  +  40*95.3 -Spc        (iv.g.i-5) 
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TABLE  VI 

KINEMATIC  AND  INERTIAL  COUPLING  FOR  CIRCULAR  AIRFRAME 

CONDTION  INERTIAL  (deg/sec2)    KINEMATIC  (deg/sec  ! 

MIN       MAX  MIN         MAX 

Circular  airframe 

with  cross-        -0.008     92.794      -1.132     0.012 

coupl ings 

(section  IV. G) 


with  cross- 
coupling  and 

faster  response     -0.0768    100.573      -2.04       0.3353 
K   =  -0.0387, 
K   =  -15.6 
(section  V.F) 


with  cross- 
coupling  and 

faster  response     -0.0212    103.553      -1.39       0.0436 
K.  =  -.0387, 
K   =  -70.88 
(section  V.F) 


witn  cross- 
coupling  and 

slower  response     -0.107     118.294      -2.97       0.0387 
K   =  -0.0274, 
KH  =  -15.6 
(section  V.F) 


with  cross- 
coupling  and 

slower  response     -0.0743    109.411      -2.32       0.0431 
K   =  -0.0274, 
KA  =  -70.88 
(section  V.F) 

with  cross- 
coupling 

slower  response    -6.0791    199.841     -0.8662     1.0541 
and  increased  gain 
in  tne  coordinated 
branch  Kyp  =  1.0, 
K.  =  -0.2§7 
K   =  -70.88 
(section  V.F) 


with  cross- 
coupling  slower 
response  and 

modified  pitch      -2.264     38.386       -0.00788    -0.6021 
control  law  as 
section  V.G 
KA  =  -0.0274 

Note:   KA,  K,  K   as  referred  in  section  V.F  and  V.G. 
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Figures  5.27  and  5.28  show  the  minimized  inertial  and 

kinematic  cross-coupling,  with  feedback  the 

angle-of-attack,  with  Ka  =  -0.0274  and  K  =  -70.88. 

Al 
Comparing  the  values  of  the  above  minimized  inertial  and 

kinematic  coupling,  with  the  values  of  the  kinematic  and 

inertial  coupling  when  the  pitch  channel  has  feedback  the 

pitch  angular  rate  (q)  and  KA  =  -0.0274,  as  given  in  table 

VI,  it  is  concluded  that  the  couplings  when  the  pitch 

channel  has  feedback  angle-of-attack  (a)  and  the 

aforementioned  particular  compensator  (Fgure  5.26),  were 

reduced  228.29%  for  inertial  and  190%  for  kinematic. 

Also,  comparing  the  Figure  5.25  which  illustrates  the 

achieved  inertial  acceleration  nz  for  circular  airframe 

with  minimized  coupling  effects  and  Figure  5.29  illustrated 

n7  for  circular  airframe  with  minimized  coupling,  it  is 

concluded  that  as  the  coupling  reduces,  the  coupling 

effects  increases. 

H.    CONCLUSIONS  AND  RECOMMENDATIONS 

The  conclusions  that  follow  are  based  on  a  single 
representative  flight  condition  (Mach  3.95  and  60  kft 
altitude)  and  for  the  particular  analysis  which  was  done  in 
this  chapter. 

1.   From  the  studies  in  this  chapter,  it  is  concluded: 
a.   Transients  in  maneuver  plane  acceleration  are 
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Fig.  5.27  Minimized  inertial  Coupling  (P.r)  vs  Time  (t) 
CBTT  of  Circular  airframe 
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Fig.    5.28      Minimized    Kinematic    Coupling    (-P.B)    vs    Time    (t) 


CBTT  of   Circular    airframe 
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Fig.  5.29  Achieved  Inertial  acceleration  (rO 
Modified  pitch  channel  for  minimiza 
coupling;  CBTT  of  Circular  airframe 


vs 
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caused  by  kinematic  and  inertial  coupling  between  pitch  and 
yaw  dynamics  through  missile  roll  rate.   Transients  are  in 
the  form  of  overshoots  and  undershoots  (section  V.B)  . 

b.  Using  as  feedback  in  the  nonlinear  pitch 
channel  for  both  airframes  the  pitch  angular  rate  (q) ,  the 
transients  in  the  maneuver  plane  exceeds  the  denoted 
requirements  and  the  achieved  maneuver  acceleration  is  not 
acceptable  (sections  V.D  and  V.F). 

c.  Using  as  feedback  in  the  nonlinear  pitch 
channel  for  both  airframe,  the  angle-of-attack  (a)  and 
modifying  suitably  the  pitch  control  law  (Figure  5.5  and 
Figure  5.20) ,  the  transients  are  minimized  below  the 
desired  requirement  and  the  achieved  inertial  acceleration 
(nz)  is  acceptable.   (sections  V.D  and  V.F). 

d.  Using  as  feedback  the  rate  of  angle  of  attack 
(a),  the  whole  system  becomes  uncontrolable  because  the 
pitch  tail  incidence  (6  )  exceeds  the  limits.   (sections 
V.D,  V.F) . 

e.  Using  as  feedback  the  rate  of  angle  of  attack 
(a)  and  modifying  suitably  the  pitch  control  law  (Figure 
5.26),  the  kinematic  and  inertial  coupling  are  reduced 
dramatically,  comparing  with  their  values  when  the  used 
feedback  is  the  pitch  angular  rate,  (sections  V.D,  and  V.G) 

f.  Using  as  feedback  the  rate  of  angle  of  attack 
(a),  in  the  modified  pitch  channel  (Figure  5.26)  for 
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minimization  of  kinematic  and  inertial  coupling,  the  whole 
system  becomes  uncontrolable ,  because  the  pitch  tail 
incidence  (fi  )  exceeds  the  limits,  (sections  V.E  and  V.G) 

g.   Minimizing  the  inertial  and  kinematic 
cross-coupling  effects,  the  values  of  kinematic  (-pB)  and 
inertial  (pr)  coupling  increase,  (section  V.E  and  V.G) 

h.   The  above  conclusions  show,  that  for  the 
particular  airframe,  flight  conditions  and  analysis  is  not 
possible  to  achieve  minimization  of  the  kinematic  and 
inertial  cross-coupling  effects  and  simultaneously 
minimization  of  the  cross-couplings  themselves. 
2.   Recommendations 

Considering  the  studies  of  this  work  and  the 
conclusions,  it  is  recommended: 

a.  Further  analysis  of  the  minimization  of  the 
kinematic  and  inertial  cross-coupling  effects  and  the 
cross-coupling  themselves  in  other  airframes,  and  with 
other  analysis  to  achieve  simultaneous  minimization. 

b.  Further  analysis  of  the  CBTT  autopilots  in 
other  flight  conditions  for  verification  of  the  above 
results. 
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APPENDIX  A 

MISSILE  SIZING  AND  MASS  PROPERTIES 
In  order  to  provide  a  realistic  missile  based  on  the 
configuration  concepts  tested  aerodynamically  in  [Ref.  1]  , 
the  models  have  the  following  characteristics. 

Characteristic  Circular       Elliptical 

Length  (in)  168  168 

Max.  Diameter  24 

Max.  Major  axis  (in)  41.57 

Max,  Minor  axis  (in)  13.86 

Weight  (lbs)  2525  2475 

I    (slug-ft2)  40  110 


xx 

yy 


I    (slug-ft2)  804  790 


I    (slug-ft2)  810  853 

eg  distance  from  L.E  (in)   100.8  100.8 

Reference  length  d  (ft)      2  2 


2 

Reference  area  S(ft  ) 


TT 
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APPENDIX  B 

NONLINEAR  AERODYNAMIC  DATA 

The  aerodynamic  data  were  taken  from  [Ref .  2]  .   For 
reference,  the  normal  and  pitching  moment  curves  at  M  = 
3.95  have  been  produced  in  fig.  B.l  and  B.2  for  the  two 
configurations  in  [Ref.  2].   The  aerodynamic  derivativies 
Cv,  C  ,  C  ,  with  respect  to  sideslip  angle  6,  yaw  control 
6„,  and  roll  control  <sD  are  presented  in  figures  B.3,  B.4 
and  B.5  of  [Ref.  2]  as  used  in  the  computer  simulation, 
namely,  as  piece  wise  linear  segments  for  ease  in 
interpolation. 
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APPENDIX  C 

REQUIREMENTS  FOR  UNCOUPLED  AUTOPILOTS 

The  requirements  for  the  classical  design  technique  of 
the  uncoupled  autopilots  are  given  in  [Ref.  2]  and  are  the 

following : 

1.  Acceleration  time  response 

a.  63  percent  time  constant  of  0.5  seconds  for  a 
step  command  of  acceleration  at  the  flight  condition  of 
interest,  (i.e  M  =  3.95  and  altitude  60kft)  and  small 
angles-of-attack.   This  response  is  representative  of  a 
tactical  missile  of  this  size. 

b.  Overshoot    10  percent 

c.  Zero  steady  state  error  in  acceleration  to 
reduce  variations  of  guidance  navigation  gain. 

2.  Relative  stability 

Gain  margins   6db,  phase  margins   30  deg  with  a 
goal  of  12  db  and  50  deg. 

3.  High  Frequency  Attenuation  in  Actuator  Command 
Branch 

It  must  be  >,   15  db  at  100  rad/sec  and  zero 
angle-of-attack.   This  requirement  will  provide  sufficient 

high  frequency  attenuation  for   30  Hz  actuator  and  for  body 
bending  modes,  when  high  frequency  filters  are  added.   This 

requirement  also  limits  autopilot  speeds  of  response. 
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